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Physicists become acquainted with special functions early in their studies. Consider 
our perennial model, the harmonic oscillator, for which we need Hermite functions, 
or the Laguerre functions in quantum mechanics. Here we choose a particular number 
theoretical function, the Riemann zeta function and examine its influence in the realm of 
physics and also how physics may be suggestive for the resolution of one of mathematics' 
most famous unconfirmed conjectures, the Riemann Hypothesis. Does physics hold an 
essential key to the solution for this more than hundred-year-old problem? In this 
work we examine numerous models from different branches of physics, from classical 
mechanics to statistical physics, where this function plays an integral role. We also 
see how this function is related to quantum chaos and how its pole-structure encodes 
when particles can undergo Bose-Einstein condensation at low temperature. Throughout 
these examinations we highlight how physics can perhaps shed light on the Riemann 
Hypothesis. Naturally, our aim could not be to be comprehensive, rather we focus on 
the major models and aim to give an informed starting point for the interested Reader. 
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I. INTRODUCTION 



'Can you do Addition?' the White Queen asked. 
'What's one and one and one and one and one and 
one and one and one and one and one?' 
'I don't know, ' said Alice. 'I lost count. ' 
(Lewis Carroll - Through the Looking Glass) 

Counting, in the broadest sense, is probably the oldest 
mathematical activity and not even uniquely ours. Even 
animals can distinguish one, two and three, maybe just 
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by recognising a pattern, but only humans have devel- 
oped an abstract language, mathematics or more specifi- 
cally number theory, which accurately describes the prop- 
erties of numbers. 

In the following we will focus on the border between 
physics and number theory, and more precisely, how the 
Riemann-zeta function, £(s), appears in quite different 
areas of physics. This review does not intend to be com- 
prehensive, rather would like to offer a panoramic view 
and give a feeling as to why many physicist find beauty 
in the structure of this seemingly random function and 
what one might learn from it. We collect examples from 
diverse realms of physics, from classical mechanics to con- 
densed matter physics, where the Riemann-zeta function 
or its 'descendants' play a significant role. Due to space 
limitations we do not aspire to be mathematically pre- 
cise in our derivations, but we give physical arguments 
to support results, and also direct the Reader to relevant 



II. HISTORICAL BACKGROUND AND 
'MATHEMATICAL NECESSITIES' 

God invented the integers; 
all else is the work of man. 
(Leopold Kronecker) 

Natural numbers form the basis of our arithmetic, with 
various operations defined among these numbers. All of 
us learn to use four basic operations: addition, subtrac- 
tion, multiplication and division. The latter, division, 
hides one of the most enigmatic internal structures of the 
set of the natural numbers, namely that there are special 
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numbers, the primes, among the natural numbers which 
cannot be divided by any other natural number, other 
than unity and themselves, without a remainder. Euclid 
of Alexandria proved that there are infinitely many such 
numbers. Later, Eratosthenes of Cyrene gave a theoret- 
ical algorithm, a sieve, for finding these primes amongst 
the natural numbers. Despite all efforts in the last two 
thousand years, the efficient determination as to whether 
a given number is prime or not still proves a remarkable 
challenge. 

It is not hard to understand why the distribution of 
primes could captivate the imagination of many mathe- 
maticians and physicists. These numbers seem to obey 
two contradictory principles. Firstly, they seem to ap- 
pear randomly among composite numbers, but secondly 
they also appear to obey strict rules governing their dis- 
tribution. 

Apart from Euclid's, numer ous proofs exist fo r the in- 
finitude of the prime numbers ( Ribenboim . 199lh . Euler, 



at th e early age of 30, proved a stronger statement ((Euler 
1737h . 
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This formula clearly proves Euclid's statement but it also 
demonstrates the frequent occurrence of prime numbers 
amongst composite numbers. A natural continuation of 
his work was the analysis of the arithmetic properties 
of the series, ^n _fc . Substituting k = 1 into this ex- 
pression we recover the well-known, divergent harmonic 
series. Conversely, if k > 1 t he summation converges. 
Euler also showed (|EulerL 117371 ) - using the fundamental 
theorem of arithmetic - that this series can be written as 
an infinite product over the prime numbers, p, such that 



(2) 



One may interpret through this relationship that the 
prime numbers construct the function. Since p de- 
notes a prime number and k > 1 , none of the factors in 
this product can be zero. Therefore we can conclude that 
does not have any zeros if k > 1. 
Bernhard Riemann, who was the first to apply the tools 
of complex analysis to this function, pro ved that the func - 
tion defined by the infinite summation ( Riemann . 18591) 
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can be analytically continued over the complex s plane, 
except for s = 1. This analytic continuation of the func- 
tion is called the Riemann-zeta function. Here we follow 
the traditional notation, with s denoting a complex num- 
ber, s — a + it, where a and t are real numbers and i is 
the usual imaginary unit. 



FIG. 1 The 'anatomy' of the Riemann-zeta function on the 
complex s plane. The black dots (•) represent the zeros of 
C(s), including possible zeros which do not lie on the critical 
line. 
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Riemann also derived a functional equation, containing 
the £(s) function, which is valid for all complex s and 
exhibits mirror symmetry around the a — 1/2 vertical 
line, called the critical line, such that 
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C(l - .) (4) 



One should note that the zeta function stands on both 
sides, on the left hand side with argument s, while on the 
right hand side with (1 — s). This relationship between 
£(s) and C(l — s) provides some insight regarding the 
location of the zeros of this function. Let us examine the 
half-line for which a < 0, and t — 0. The products on 
either side can be zero if at least one of the factors is zero. 
On the right hand side of (@J all the pre-factors of the zeta 
function are non-negative and do not have any zeros. On 
the other side, however, the T(a/2) function has simple 
poles at all even negative integers. The equation can 
hold only if ((a) has simple zeros at the same locations. 
These zeros are called trivial, because their locations are 
inherited from the T function. The same argument also 
shows that all other zeros of the £(s) function have to lie 
in the < a < 1 region, called the critical strip. The 
zeros located in this strip are the non-trivial zeros of the 
Riemann-zeta function. It can also be shown that the 
non-trivial zeros p are arranged symmetrically, both in 
respect of the critical line and the t = axis. Figure Q] 
depicts the pole and zero structure of £(s) on the complex 
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s plane including the possible zeros off the critical line. 

So far the statements about the zeros of £(s) and their 
locations on the complex plain were simple. However 
the distribution of the non-trivial zeros holds one of the 
most intriguing and enigmatic mathematical mysteries of 
the last century and a half. It is embarrassingly easy to 
pose Riemann's conjecture: all non-trivial zeros of £(s) 
have the form p = 1/2 + it, where t is a real number. In 
other words all non-trivial zeros lie on the critical line. In 
1900 Hilbcrt nominated the Riemann Hypothesis as the 
eighth problem on h is famous list of compelling prob- 
lems in mathematics (jHilbertl . Il902l) . Since then not just 
professional mathematicians but mathematical soldiers 
of fortune tried, and still try, to verify its validity. The 
stakes are high. Whoever proves or disproves this hy- 
pothesis engraves his name in the tablets of the history 
of mathematics, and may also receive one million dollars 
from the Clay Mathematics Institut^- 

During the past century, the Riemann Hypothesis has 
been recast into many equivalent mathematical state- 
ments. A few of them are purely number theoretical 
in origin, such as the Mertens conjecture, which we will 
later discuss in the context of a special Brownian motion, 
but other redefinitions are very much cross-disciplinary. 
A more advanced mathematical introduction to the his- 
tory of the Riemann Hypothesis and its equivalent state- 
ments ca n be found in an exce llent monograph and com- 
pendium ( Borwein et all [20081) which is readable not just 
at the expert, but also the undergraduate level. 

The distribution of the £(s) zeros, with real part equal 
to 1/2, has thus attracted significant interest. One of 
mathematics' giants has pro ven tha t infin itely many ze- 
ros do lie on the critical line ( Hardvl . Il9l4) , however Rie- 
mann's conjecture is much stronger, requiring all the ze- 
ros to be on the critical line. In 1942 Selberg proved 

N (T) > C Tln(T) (C > and T>T ) (5) 

i.e. the number of zeros of the forms = \+it (0 < t < T), 
denoted by Nq(T), grows as Tin (T) at least for large T. 
Three decades later, in 1974, Levinson showed that at 
least one third of the non-trivial zeros are on the critical 



line (jLevinsonl . 119741) which was la ter incrementally im- 
proved to two fifths ( Conrev , 1989h . This small step over 
a period of twenty years is indicative of the difficulty of 
the Riemann Hypothesis. 

Let us return to the linkage between the £(s) zeros 
and prime numbers. Equation ([2]) clearly shows the 
strong connection between the £(s) function and the 
prime numbers. This relationship can be made even more 
explicit if one examines how the number of primes be- 
low a given threshold behaves as this threshold is in- 
creased. Based on empirical evidence, man y mathe- 
maticians, e.g. Legendre, Gauss, Chebyshev (jDicksonl . 



20051) . have conjectured that the prime counting func- 
tion, n(x) = | {p | pis prime and p < x}|, asymptotically 
behaves as the logarithmic integral Li(x). This conjec- 
ture is known nowa days as the Prim e Number Theo- 
rem afte r Hadamard (iHadamardl 18961) and de la Vallee- 
Poussin ( de la Vallee-Poussinl 18961) independently gave 
rigorous proofs of this statement. Interestingly, this theo- 
rem has a geometrical interpretation: the Prime Number 
Theorem is equivalent to the assertion that no zeros of 
£(s) lie on the a = 1 b oundary o f the c ritical strip. 

Riemann published ( Riemann . 18591). although Man - 
goldt provided the rigorous proof (jvon Mangolddll895h . 
the following explicit formula for the prime-counting 
function 7r(x) 



(6) 



where 



J(x) =Li(x) 



lim 

T->oc 



J2 EH(plog(aO) 

\P\<T 



dt 



log (2). 



(*»-_.)* log (i) 

Here p(n) is the Mobius function^, p denotes the non- 
trivial zeros of the Riemann C(s) function, and Li(x) 
and Ei(x) stand for the logarithmic and exponential in- 
tegral^, respectively. Therefore, whoever knows the dis- 
tribution of the non-trivial zeros of £( s ), wm a l so know 
the distribution of the prime numbers. 

Selecting only the first terms of the summand in equa- 
tion ([5]) reproduces exactly the Prime Number Theorem, 



7r(x) = Li(x) 



ln(x) 



(7) 



This observation may lead us to conclude that Li(x) gives 
the main contribution to 7r(x) while the other terms rep- 
resent corrections, similar to a perturbative calculation 
in physics - an analogy to which we will return. Figure 
[2] depicts the prime counting function, 7r(x) and its var- 
ious approximations. One may notice that the leading 



2 The Mobius function is defined as follows: fJ.(l) = 1, p.{n) = if 
n has a square divisor, and fj,(p±p2 ■ ■ ■ Pk) = ( — l) fe if all Pi s ar ^ 
different. Thus fi(2) = -1 and ^(12) = 0, and (1,(21) = 1. 

3 The notation for the logarithmic integral is ambiguous in the 



literature. There arc two definitions 



h(x) 



Jo 



dt 
ln(x) 



and l2(%) 



i: 



dt 
InO) 



1 See http://www.claymath.org/millennium/Riemann_Hypothesis 



where I\ is interpreted as a Cauchy principal value. These inte- 
grals differ only by a constant number. Depending on the book 
the Reader may consult, either Ii(x) or l2(%) is denoted with 
hi(x). Here, we prefer the former. 
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term, Li(x), captures the tendency of ir(x) well and the 
appearance of the oscillations clearly show how the zeros 
p n influence and refine the agreement. As x — > oo the 
curves of Li(a?) and ~k{x) will practically coincide on a 
similar plot. 




FIG. 2 Figure depicts the approximation ((6| for the prime 
counting function, n(x) (dashed line), using only the first 
term, Li(:r) (dash-dotted line), and using the first ten non- 
trivial pairs of zeros of the Riemann £(s) function (solid line). 
In the inset we restricted the range to [2, 10] and used the 
first 30 non-trivial zeros. 

One may define a density for the complex, non-trivial 
Ricmann-zeta zeros as 



d(N) =J2 S ( N ~ Pk) 



(8) 



where S is the Dirac-delta distribution. Following Sir 



Michael Berry (jBerrvl 119851 ) the spectral density can be 
separated into a smooth and an oscillatory part, d(N) = 
d{T) + d osc {T), as 



1 ( T 
d(T) = — In — 



d osc (T) 



_^ OO 

EE 



+ 1-- + 0(T-) 
In (p) cos (rT In (p) ) 



r=l 
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(9a) 



(9b) 



where the external summation of d osc (T) runs over the 
prime numbers, p. The oscillatory part, therefore, gives 
the fluctuations as individual contributions from each 
prime number p labeled by an integer r corresponding 
to the prime power p r . Based on the smooth density of 
Riemann-zeros one may derive the number of positive, 
non-trivial zeros upto a fixed value of Tq: 



N(p < T ) 



T 



d(T) dT 



2ir 



In 



To 
2tt 



(10) 



Changing variable to T ■ 
result using T we obtain 



ln(To/27r) and recasting our 



i.e. the number of £(s) zeros below T increases expo- 
nentially. Although at this point this change of variable 
seems somewhat arbitrary, we will see later that it fur- 
ther strengthens the similarity between the zeros of £(s) 
and the periodic orbits of a chaotic system, where the 
number of periodic orbits also increases exponentially. 

Finally, we note the fruitful and diverse area of exten- 
sions of the Riemann-zeta function. These generalised 
zeta-functions do also occur throughout physics, primar- 
ily in modern quantum field theories. This topic, how- 
ever, is far beyond the scope of this shor t review a nd we 
can only suggest Elizalde's monograph (Elizaldc, 1995) 
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as an introduction and Lapidus' book ([Lapidus 
for a more authoritative study. 



III. CONNECTIONS TO PHYSICS 

The Riemann Hypothesis is a precise statement, 
and in one sense what it means is clear, but 
what it's connected with, what it implies, where 
it comes from, can be very unobvious. 



(Martin Huxley) 



A. Classical mechanics 

In this section we discuss those models of classical me- 
chanics, such as billiards, which lead to the introduction 
of the notion of integrability and chaos. This develop- 
ment of ideas gave birth to a new paradigm, since it 
provided an insight into how the spectrum of quantised 
analogues of classical systems are connected to classical 
paths. 

Classical mechanics, in its Lagrangian and Hamilto- 
nian forms, is the exemplar for physics in the modern 
sense. The major theories, e.g. statistical mechanics, 
quantum mechanics, are first expressed in the language 
of analytical mechanics with the development traced to 
the Enlightmcnt. Although a few analytically solvable 
models, e.g. Kepler two-body problem, harmonic oscilla- 
tor, gave confidence in the machinery of mechanics, it was 
soon realised that there are important cases, e.g. three- 
body problem, where one not just cannot solve the equa- 
tions of mo tion analytically, but the m otion is proven to 
be chaotic ( Celletti and Perozzii l20C)7h . This behaviour 
is very peculiar and at first sight seems puzzling, since the 
governing equations are deterministic, yet the actual mo- 
tion seems to behave randomly. The celestial relevance 
of this three-body problem was so fundamental and en- 
ticing that King Oscar II of Sweden and Norway offered 
a pri ze for the person who could solve the following prob- 
lem ( Barrow-Greenl 11994 ) 



N{T) oc e 



r 



(11) 



For an arbitrary system of mass points which 
attract each other according to Newton's law, 
assuming that no two points ever collide, give 
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the co-ordinates of the individual points for 
all time as a sum of a uniformly convergent se- 
ries whose terms are made up of known func- 
tions. 

Although this problem had not been solved, Poincare was 
awarded this illustrious prize for his impressive contribu- 
tion. His work revolutionised the analysis of such chaoti- 
cally behaving systems, although one had to wait nearly 
a hundred years for this revolution to really happen. 

In classical mechanics we distinguish a special class 
of systems, the integrable dynamical system, which pos- 
sess as many independent integrals of motion, X n , (action 
variables) as degrees of freedom, N. For these systems 
the Hamiltonian can be expressed as a function of these 
action variables, namely % = Will, . . . ,1^), and the 
equations of motion (n = 0, 1, . . . , N) 



dcp 



dH , dl n 
— 7t~t~ and — — 

at oln dt 



dH 



(12) 



are easy to solve: I n — constant and cp n — ip n fl + co n t. 
A theorem of topology then guarantees that these TV 
constants of motion, provided they are independent of 
each other, define an N dimensional torus and each tra- 
jectory with constant energy lies on that torus. There- 
fore, as a specific case, the dynamics described by a one- 
dimensional time-independent Hamiltonian is necessarily 
integrable. In order to consider chaotic dynamics one has 
to either introduce a time-dependent Hamiltonian or in- 
crease the degrees of freedom to two or higher. 

One of the 'simplest' generic models with two or more 
degrees of freedom is that of classical billiards. These 
are dynamical systems where a particle has constant en- 
ergy and moves in a finite volume, which may contain 
impenetrable obstacles. Whenever the particle reaches 
the boundary it suffers specular reflection. Depending 
on the shape of the billiard, the motion can be integrable 
or chaotic. The analysis of a circular billiard (see Fig- 
ure [3]) is straightforward due to the rotational symmetry 
The incident angle remains the same at each bounce and 
each impact can be calculated from the previous one by 
rotating the circle twice that angle. Therefore if the in- 
cident angle is a rational multiple of n, i.e. mn/n, the 
trajectory is periodic with period n and therefore finite, 
otherwise it is infinite. In this latter case the points where 
the ball hits the wall will be uniformly distributed along 
the circumference of the circle. It was also proven by 
Jacobi that in the latter case every interval of the circle 
contains points of the trajectory. 

Before we step beyond billiards and generalise the 
idea of periodic orbits, the origin of trace formulae, 
let us make a short detour aro und a recent result 
( Bunimovich and Dettmannl . 12005 ) regarding the circu- 
lar billiard (see Figure HJ. As we discussed, due to 
rotational symmetry, or in other words, the conserva- 
tion of angular momentum, this billiard model is inte- 





FIG. 3 A circular billiard and a Bunimovich stadium, which 
is a rectangle smoothly joined by semi-circles. Two differ- 
ent types of trajectories, periodic orbits (1) and non-periodic 
trajectories (2) are also depicted. 



grable and the trajectory is fully described by two an- 
gles, P and ip, the angle around the circumference mea- 
sured from a pre-determined point and the incident an- 
gle of the trajectory at the boundary, respectively. With 
these variables the dynamics is governed by the map- 
ping: (/3, n- ((3 + n — 2ip,ip), where all angles are 
taken modulo 2ir and the ball travels with unit veloc- 
ity. The phase space of this system can be described 
by Birkhoff's coordinates constructed from two angles; 
the arc-length coordinate q = f3 (measured in radians 
and modulo 2tt), and the tangential momentum coordi- 
nate defined as p — sin(^) . By convenient normalisation, 
the arc-length of the billiard is unity and the velocity 
of the ball is also unity, the phase space is restricted to 
< q < 2-7T, and —1 < p < 1. This choice also introduces 
a natural unit time-step, the time elapsed between con- 
secutive bounces, At = 2 cos (?/>). The movement of the 
ball can, therefore, be represented by a possibly infinite 
series of points inside this phase-space area. Despite the 
rather artificial appearance of this model, the electromag- 
netic field in opti cal or microwave cavities can be modeled 
by such billiards (Alt et al\ 1998 : Haravama et al , 2003 ; 
Nockel et ail Il997t IStockmann and Steinl . Il990lh Since 
these experimental billiards are not ideal, it is interest- 
ing to examine what happens to the dynamics of this 
system if we cut a small window(s) along the reflective 
boundary, thereby, naturally introducing dissipation or 
'leakage'. It is natural to ask: what is the probability, 
P(n), of a ball leaving the billiard after n bounces, what 
is the mean number of bounces, (n), before the ball es- 
capes, or similarly, what is the probability, P(i), that 
escape takes at least time t. 

For strongly chaotic billiards the latter probability de- 
cays exponentially, while for integrable billiards, such 
as the circular one, it soft ens to only power-law decay 
((Bauer and Bertsch I . |l990) and can be qualitatively un- 
derstood using a simple geometrical argument. The prob- 
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ability, p, that the ball escapes in a bounce is proportional 
to the size of the gap to that of the boundary, p = e/L, 
Moreover, the probability that the ball survives the first 
(n. — 1) bounces and escapes only at the nth bounce is 
(1 - p) {n -^p. Therefore the mean number of bounces 
occurring until escape is 



\ncscapc, 



= £fc(i-p)^p=i«i 
fc=i y 



(13) 



all the integer numbers smaller than n, one could con- 
nect it to the Fourier series. However, here one only uses 
those m's which are relative primes to n. Converting the 
cosine term to complex exponentials and using Ramanu- 
jan's identitjd for the sum of exponentials, the contribu- 
tion of the cosine term turns out to be another special 
function of number theory which we have already met, 
the Mobius function, fJ,(n). Therefore, the probability of 
non-escaping orbits is 



Let us now cut two (possibly overlapping) holes, with 
sizes e, on the boundary and examine the non-escaping 
periodic orbits. Based on the geometrical argument used 
above, we expect the probability to be ~ 2/e, if the two 
holes do not overlap. However, in systems where the tra- 
jectories do not diverge strongly, i.e. Lyapunov exponent 
is close to zero, only a small fraction of the trajectories 
will eventually hit the opening on the boundary, and the 
mean escape time will be proportional to e. 

If the initial incident angle is taken to be ipm,n = 7r/2 — 
rmr/n, where m < n are integers and relative primes to 
each other, then the trajectory is closed and its period 
is n. Let us now examine only those initial conditions 
for which the escape time is at least t, or in other words, 
the number of bounces is at least N = \ 2ir/e\ . To fulfill 
this requirement one might take the initial value of ip = 
ipm, n + where < rj -C e and /3 can be restricted to 
the following range 



nt 



COS (lpm,n) 



u 



Tft 



2n 



COS (lp m ,n) n 



The prime indicates that angles are taken modulo 2ir/n. 
The probability can, therefore, be calculated if one sums 
up all possible values of (m, n) pairs. This is the point 
where number theory enters into this physical problem; 
we have to guarantee that m and n are relative primes. 
Integrating over the permitted region of /3' one may find 



P(i,e,0)~i$>.nn)£ 



1 — cos 



/ 2rmr 
V n 



(14) 



where the exact form of ^F(n) can be found in 
( Bunimovich and Dettmannl . l2005f) . Surprisingly the 
sum over m can be explicitly determined. The first, 
unit term, simply counts how many numbers are relative 
prime to n and, therefore, it can be formally expressed 
using a special function of number theory; Eulcr's to- 
tient function 0. The second term in equation (|14j) is 
also a special expression. If the summation were over 



4 Euler's totient function, (j>(n), gives the number of positive inte- 
gers smaller than n, which are relative prime to n, e.g. for any 
prime number <p(p) = p — 1, since all integers smaller than p are 
relative prime to p. 



Poo = lim (tP(t,e,0)) 



n=l 



n [cj)(n) - n(n)} T{n) (15) 



The leading order behaviour of Poo as a function of e can 
be determined by calculating its Mellin-transform 



P(s) = / P 00 (e,e)e s - 1 de 



(16) 



and examining the residues of P(s) on the complex s- 
plane. Bunimovich and Detteman showed that for the 
two-hole problem, where these holes are separated by 0°, 
60°, 90°, 120°, 180°, the probability P(s) is uniquely 
determined by the Riemann-zeta function, C( s )i by i.e. 
its pole and non-trivial zeros. The first corrections to 
the leading order term are given by the non-trivial zeros 
of C(l + s), which are of the order y/e ln(e) m_ , provided 
for all zeros of C( s )> 3?(s) = a < i with multiplicity m. 
The Riemann Hypothesis is then shown to be equivalent 
to different asymptotic estimates on t he number of zeros 
(|Titchmarsh and Heath-Brownl . 120031) . Therefore, if the 
non-trivial zeros provide the second order corrections to 
the probability, it is instructive to examine the deviation 
of these probabilities experimentally from the leading- 
order geometric terms, namely 



lim lim ( e 5 - 1/2 



>ot- 



tPi(t) - - 
e 



lim lim [ e 5 - 1 ' 2 



tPi(t)-2tP 2 (t) 



(17a) 



= (17b) 



where Pi, P2 belong to the one- and two-hole problem, 
respectively. If it is (experimentally) found that for ev- 
ery 8 > these equations are fulfilled, then it proves 
the validity of the asymptotic formulae, thus the valid- 
ity of the Riemann Hypothesis. The numerical results 
by Bunimovich and Detteman does not contradict these 



5 Ramanujan's sum is defined as 
Cn(m) 



where the summation is over those values of m, which are relative 
prime to n. Using Mobius inversion fo r this sum one can prove 
that c n (m) = ,u(n) (Hard y and Wrighdll960h . 
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It can b e shown rigorously for a generic classical chaotic 
system ( Cvitanovic and Eckhardti Il99lh . that 



FIG. 4 A circular billiard w ith to small openings 
(jBunimovich and Dettmannl . I2005T ) . 



equations. Although their result has not proven the Rie- 
mann Hypothesis, it provides a physically realisable sys- 
tem where actual measurements can substantiate, but 
not prove, the conjecture. 

Let us turn our attention now to the dynamics 
of a more general billiard system. If one smoothly 
deforms the boundary of this circle and creates a 
stadiu m-like shape , the a nalysis is far less straightfor- 
ward ( BunimovichL Il979t ). However, qualitatively we 
may see that the trajectories can be classified sim- 
ilarly and one may distinguish periodic orbits, i.e. 
{q n {h),Pn(ti)} = Un{t2),Pn{h)} for ti < t 2 , and non- 
periodic trajectories. It is tempting to think that periodic 
orbits are exceptional and quite rare among all orbits, 
since for integrable systems the number of periodic or- 
bits grows polynomially and one may expect that violat- 
ing integrability would decrease the number of periodic 
orbits. In fact, the opposite is true. These special or- 
bits proliferate among the possible orbits and their num- 
ber, for a general Hamiltonian dynamics, grows exponen- 
tially with the length of the periodic orbits, ~ e /hi, 



where h is called the topological entropy an d I denotes 
the length of a given cla ss of periodic orbits ( Gutzwillen . 
199lt IStockmannl Il999t) . This is a striking difference be- 
tween integrable and chaotic systems. It is even more sur- 
prising that the knowledge of these periodic orbits serves 
as a powerful analytical tool for investigating chaotic sys- 
tems, and, moreover, they provide the pathway, through 
trace formulae, from classical to quantum mechanics. 

Let us examine the time evolution of a Hamiltonian 
flow in general. We denote the trajectory starting its time 
evolution from the initial point ro with r(t) = F(ro,£). 
We further introduce the evolution operator with the fol- 
lowing definition 



Tr(£(i;r',r)) =^T p ]T 



S(t - rT p ) 

rj |det (i-J;) 



(19) 



C{t- r', r) = <5(r' - r(t)) = 6{r' - F(r, t)) 



(18) 



where the first summation runs over the periodic orbits 
labeled by p, while the second takes into account all repe- 
titions, r. J p is the Jacobian matrix of F localised around 
the periodic orbit, also called the monodromy matrix. 

Here we can make an important observation: although 
this equation looks cumbersome, it does relate the spec- 
trum of the evolution operator to a global behaviour of 
periodic orbits. Therefore these two sets of abstract ob- 
jects are intimately related to one another. The connec- 
tion of this trace formula and its quantum mechanical 
counterpart to the Riemann zeta function will become 
clear in the next section. 



B. Quantum mechanics 

Below we expound the Polya-Hilbert conjecture. We enu- 
merate the one- dimensional Hamiltonians proposed for 
which the distribution of energy eigenvalues mimic the 
non-trivial zeros of the Riemann zeta function and anal- 
yse their relationship with the Gutzwiller trace formula. 
We also examine the possible symmetries of a 'Riemann- 
operator' since it partially encouraged the development of 
quantum mechanics with only CT or PT symmetry. 

In the dawn of the 20th century Bohr postulated a se- 
ries of rules for describing the spectrum of the hydrogen 
atom well before the birth of Schrodinger's and Heisen- 
berg's quantum mechanics. In these early days 'quantisa- 
tion' meant to restrict the possible values of action vari- 
ables of the classical system (Bohr-Sommerfeld, Wentzel- 
Kramers-Brillouin, etc.) and the rules worked well, upto 
an additive constant. However, this description cannot 
be satisfactory in general, since for the majority of clas- 
sical systems the only constant of motion is the energy, 
and therefore a method of quantisation relying on the 
existence of act ion-angle variables could not be applied 
(|Einsteinl . 1 19 171) . On the other hand, we know that clas- 
sical mechanics works well for large systems, therefore 
quantum mechanics must give the same predictions for 
a large system as classical mechanics (Bohr's correspon- 
dence principle). This unproven principle ties these two 
theories firmly together and the same principle inspired 
the use of the Riemann £ function in investigating the 
relationship of classical to quantum mechanics. 

The basic question is: how can we quantise a classi- 
cal mechanical system? Could we state anything about 
the spectrum of a quantum system, at least qualitatively, 
without solving the corresponding Schrodinger equation? 

We cannot expect to be able to infer the complete spec- 
trum of a generic system, but asking only for the aver- 
age density of states may prove feasible. One can give a 
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crude, although remarkably precise, estimate: each quan- 
tum state occupies approximately tJ phase-space, where 
H is the Planck constant divided by 2tt and / is the num- 
ber of degrees of freedom. This result is rather general 
and the individual quantum systems differ only in the 
'fluctuations' around this average. It turns out that the 
type of fluctuation depends on the behaviour of the clas- 
sical counterpart; classically regular and chaotic systems 
are quite different. For example, a classical rectangular 
billiard is intcgrable (regular), whereas its quantum ana- 
logue exhibits a chaotic spectrum, with the spacing of 
the quantum levels, s = e n — e n _i, following an exponen- 
tial distribution P(s) ~ e~ s (see later in section BlI.CI) . 
However, a stadium billiard is classically chaotic, but the 
spectrum of its quantum counterpart is regular, mean- 
ing that P(s) is small for small values of s and sharply 
peaked at a finite value indicating a regular distribution 
of energy levels. We will see that interpreting the £(s) 
zeros as energy levels their distribution is breathtakingly 
similar to those of a quantum system's. This has inspired 
physicists to examine whether one could associate a dy- 
namical system with the Riemann zeta function. 

The advantage of this approach would be that the huge 
number of £(s) zeros are known and quick numerical al- 
gorithms have also been developed to find further ze- 
ros, thus solving the Schrodinger equation for large ener- 
gies would be unnecessary. The Riemann zeta function 
could play the same role in the examination of chaotic 
quantum systems as the harmonic oscillator does for in- 
tegrable quantum systems. This is the point where the 
examination of the Riemann zeta function may help to 
understand physics or, vice versa, the physics may lead 
us to the solution of this so far intractable mathematical 
problem. 

In order to establish a strong formal connection be- 
tween a generic chaotic quantum system and the distri- 
bution of the Riemann £(s) zeros, we have to elucidate 
a new description of quantum systems, the Gutzwillcr's 
trace formula. This trace formula is the analogue of equa- 
tions (|9alb) for physical systems. 

Let us, therefore, return to a classically integrable sys- 
tem, for which the Hamiltonian % can be given in terms 
of conserved quantities % = H(L\,... ,2jv)- Using Bohr's 
semiclassical quantisation rules, these action variables 
take not arbitrary, but fixed values 



semi-classical approximation 



Ik. 



\n k + - 



(fc = l,2,...,iV) 



(20) 



where the Uk are integers and called Maslov in- 
dices ([Ampl'dl Il997l ). The density of states, therefore, 
becomes 



d(E) = J2s(E-n(i)). 



(21) 



which can be recast as the sum of a smooth and an oscilla- 
tory term. The former originates from the Thomas-Fermi 



d TF (E) = J 6(E-H(jp,q))^ 



dpdq 



(22) 



while the latter is obtained by expanding the effective 
action to quadratic order around the classical p eriodic 
orbits (jBerrv and Taborl . Il976t lEmile et all 120061 ): 



N 



(/-3)/2 



ftV det (NQijN) 



x exp - i-Nfx + i-/3) , (23) 
\ n 4 4 / 

where Qjj = det (H) x %r. is the co-matrix of Hi,j = 
difii-H, while (3 is related to the signature oiHij. 

We see, as in classical mechanics, one can also ex- 
press the density of states as a sum of a smooth function 
dTF{E) and an oscillatory function which is defined on 
the periodic orbits of the semiclassical system. Due to the 
correspondence principle, we expect the Thomas-Fermi 
density of states to remain valid and only the oscillatory 
part to vary compared to the semi-classical derivation. 

For non-integrable systems, however, the orbits no 
longer lie on invariant tori and a different method is 
needed for the evaluation of the trace 



d(£) = --Tr(3(G B (r,r))) 

7T 



(24) 



where G.e(r, r) is the Green- function associated with a 
given Hamiltonian %. This new ap proach, based on 
the G reen-function, was developed by (jGutzwiller , Il970l 
19711 ). Here we shall not follow the details of the deriva- 
tion, but only present the final, fully quantum mechanical 
expression for the density of states 



dosc(E) 



rp 00 COS i 
± P \ ~* 



ft 2^P 



pTo. 7 ^^! |det(M«-l) 



1/2 



(25) 



where the summation runs over all primitive periodic or- 
bits, and M p is the monodromy matrix for these prim- 
itive periodic orbits. Using this new method one can 
derive a semiclassical expression for the spectrum of 
a quantum system whose classical analogue is chaotic, 
when the usual Bohr-Sommerfeld quantisation rules can- 
not be applied. Gutzwiller's result above, therefore, can 
be viewed as a bridge between the classical and quan- 
tum behaviour of a system, and can provide a rule as 
to how to quantise such a system. In this interpreta- 
tion, Gutzwiller's approach is similar to Feynman's path 
integral description, where the quantum system is de- 
scribed in terms of an infinite sum over classical paths. 
For the interested reader we can suggest, without any 
reservation, Gutzwiller' s comprehensive book on classi- 
cal and quantum chaos (jGutzwillerl . 119911) and Brack and 
Bhaduri's monograph giving an overview of semiclassical 
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physics ([Brack and Bhaduril I2003T) . In order to help the 
reader to visualise the emergence of periodic orbits in 
a quantum mechanical system we reproduce here a few 
quantum 'scars' from Heller's numerical study. Figure 
[5] shows the probability distribution for three quantum 
eigenstates of the Bunimovich billiard. It is apparent how 
the isolated, unstable classical periodic orbits manifest 
themselves as paths along which the probability distri- 
bution is greatly enhanced. Gutzwiller's idea to extract 
eigenvalues of a chaotic system via the periodic orbits, 
therefore, seems most plausible. 



TABLE I Dictionary for translating the 'Rie- 
mann dyna mics' onto a cha o tic qu a ntum dynamics. 




FIG. 5 Three eigenstates of the quantum stadium billiard are 
shown together with the major contributing unstable periodic 
orbits of the classical counterpart as thick solid lines. In the 
middle figure the guiding straig ht line for th e A shaped pe- 
riodic orbit is omitted. From ijHelkd . 119841 ) with the kind 
permission of the author. 



Based on analogy between the oscillatory part of 
the semiclassical density of states (|25|) and the sim- 
ilar expression of (l9bl) one can set up a dictio- 



nary ( Berry and Keating! 1999b : Bohigas . 2005 ) which 
maps the Riemann zeta function onto a so-far unknown 
chaotic quantum mechanical system. 

Although the Hamiltonian, "H, which would describe 



qtiics onto a cna o tic qu a ntum dynamics . 
Based on (|Berrv and Keatind . Il999bl ; iBohigasl 120051 : 
iBrack and Bhaduril . I2003D . 



Generic chaotic 
system 



Riemann zeta 
function 



periodic orbit labels 
dimensionless action 
periods 

stability factor* 
Maslov index* 
asymptotic limit 



integers 

s v /h 

det (M™ 
7i-> 



primes 
Tln(p) 
In (p) 

P r 
2t 

T v -> oc 



* Depending on how one maps the oscillatory part of the zeta 
zeros density l|9b[l onto Gutzwiller's trace formula (1251) the 
definition of the stability fa ctor and the Maslov index can be 
different. Here we followed (Brack and Bhaduri, 2003), while 
another mapping can be found in teerrv and Keatind . Il999b0 . 

' Therefore the Maslov phase is 7r, but this is not unique and 
one could also choose 3n, bit, etc. 



the chaotic quantum system corresponding to the Rie- 
mann zeta function, is still missing, the mutual resem- 
blance of (l25l) and (l9bl) reveals some p ossible properties 
of H. In (jBerrv and KeatingL 1999b ) a thorough and 
concise summary of these properties can be found from 
which we cite but a few for later use: 



1. H has a classical counterpart, since the absence of 
any analogue of h from (|9b|) indicates the scaling of 
the dynamics, namely the trajectories are the same 
at all energy scale. 

2. The Riemann dynamics is chaotic and unstable. 

3. The dynamics lacks time-reversal symmetry. 

4. The dynamics is quasi one-dimensional, because for 
a generic d dimensional scaling system the number 
of energy eigenvalues increases as ~ E d while for 
C(s) the number of zeros T < N(T) ~ T In (T) < 
T 2 . Moreover, the appearance of yfjf in the de- 
nominator implies one expanding direction and no 
contracting one. 

Below we pursue the proposed dynamics related to the 
Riemann zeta function. 

In the early days of quantum mechanics Hilbert and 
Polya suggested a physical way to verify Riemann's Hy- 
pothesis: 

I spent two years in Gottingen ending around 
the begin of 1914. I tried to learn analytic 
number theory from Landau. He asked me 
one day: "You know some physics. Do you 
know a physical reason that the Riemann Hy- 
pothesis should be true." This would be the 
case, I answered, if the nontrivial zeros of 
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the Xi-functioijfl were so connected with the 
physical problem that the Riemann Hypothe- 
sis would be equivalent to the fact that all the 
eigenvalues of the physical problem are real. 

I never published this remark, but some- 
how it became known and it is still remem- 
bered. (Private letter to OdlyzkcQ.) 

The zeros of £(s) can be the spectrum of an operator, 
1Z = + \H, where H is self-adjoint. This operator 
T-L might have an interpretation as a Hamiltonian of a 
physical system and, therefore, the key to the proof of 
the Riemann Hypothesis may have been coded in physics. 
Since the first occurrence of this conjecture a number 
of models have been promoted. Below we separate the 
models depending on whether they relate the zeros to the 
positive energy spectrum, i.e. the scattering states of a 
physical system, or to the negative energy spectrum, i.e. 
to the bound states of a quantum system. 



1. Scattering state models 

Let us first consider the possibility that the Riemann 
zeta function is associated with a quantum scattering 
problem. 

A few decades after Riem ann created a ne w geometry 
with his revolutionary work ( Riemann , 18671 ) , Hadamard 
examined the geodesies, the trajectories of freely mov- 
i ng bodies, o n sur faces with negative curvature in detail 
( Hadamardl . 1898f ) and noticed the occurrence of fami- 
lies of geodesies whose cross-section exhibits a fractal- like 
structure, as we would call it nowadays. These geodesies 
diverge exponentially, thus the distance between two tra- 
jectories, S(t), however small initially, will grow exponen- 



tially, S(t) s» e 5(0), where A is a positive number, called 
the Lyapunov exponent. This sensitivity of the system to 
the initial conditions, however, would not necessarily re- 
sult in chaotic behaviour, provided the space for the tra- 
jectories is infinite. However, if the surface is compact, 
the trajectories cannot escape to infinity, rather mix on 
this surface. If one wishes to visualise a particular exam- 
ple, consider a donut with two holes. On this surface the 
trajectories remain bounded on the s urface without the 
length of a geodesies being l imited ( Balazs and Vorosl 
19861 : iBogomolnv et all Il995t iGutzwillerl . Il99ll ). These 
two properties, exponential sensitivity of the initial con- 
ditions and mixing, are the main requirements for chaotic 



motion ( Cvitanovic et al , 2010f) . The relative simplicity 



of the description of such surfaces with negative cur- 
vature, and the presence of completely chaotic classi- 
cal motion motivated several authors in the mid-1980s 



6 G. Polya refers here to the Riemann f (s) function. 

7 See the scanned pages on Odlyzko's personal website: 
http : //www. dtc .umn. edu/~odlyzko/polya/ 



(|Balazs and Vorosl Il986t iBerrvl Il987t IGutzwillerl . Il983l ) 
to examine how such a system can be quantised, i.e. what 
properties do the solutions and eigenvalues of the equa- 
tion Hcf) = \cj> possess. 

More precisely, for free motion, one seeks the solution 

of 



A</>„ = A n </>„ 



(26) 



where <p n are required to be square integrable and 
the appropriate boundary conditions are also provided. 
Over a compact domain equation (1261) has only dis- 
crete eigenvalues. On a surface with negative curva- 
ture, the non-euclidean Green's theorem shows that the 

eigenv alues must have the form \ n = i + ip n (p is 

real) (|Gel'fand and Piatezkii-Shapirol 19591) . This re- 
semblance immediately suggests a connection with the 
zeros of the Riemann ((s). It is also proven that, for a 
compact surface, the set of n's is finite, but for a non- 
compact surface, a continuous part of the spectrum can 
also appear. In the latter case the scattering (continuous 
spectrum) is non-conventional, because it is the result 
of the geometry (curvature, compactedness) and not the 
physical interaction between particles. 

In order to express the eigenvalue density the Green- 
function is needed. Interestingly, on a surface with nega- 
tive curvature the Green-function can be explicitly writ- 
ten as a sum of individual Green-functions corresponding 
to the periodic orbits. It is also a fact that, all peri- 
odic orbits are unstable and their action is S(E) = ki, 
where k is the momentum related to the energy by 
2mE/h 2 = k 2 + 1/4 and £ defines the length of a closed 
geodesic belonging to a given conjugacy class. In this ge- 
ometry the de nsity of s tates is expressed by the Selberg 
trace formula (jSelbergl . Il949l ) 



.4 



p(k) = — k tanh (kit) , 

27T ZTT 



^ oo 

EE 



[p] 



£ p cos (nk£ p ) 
^-J sinh (ni p /2) 



(27) 



where A is the area of the surface, the first summation 
runs over conjugacy classes of primitive elements p, the 
second, their repetitions. It is important to note, the 
Selberg trace formula holds exactly, in contrast to other 
trace formulae, because no semi-classical approximation 
has been applied, although its convergence property is 
similar to the Gutzwillcr form: for large k the Selberg 
and Gutzwiller trace formulae converge, since the met- 
ric is locally Euclidean and waves with short wavelength 
lose their sensitivity to the local curvature of the metric. 
In this system, the transient scattering states were ex- 
amined by Pavlov and Fadeev who related the nontrivial 
zeros of the zeta function to the complex poles of the 
scattering matrix (jPavlov and Fadeev . 1975): 



S(k) 



~2ik 



r(|-Hfc)C(l + 2zfc) 
r(|- ik)Q(l-2ik)' 



(28) 
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Despite this natural occurrence of the Riemann zeta func- 
tion and its non-trivial zeros, no further insight into 
the zeros has been gained via this route. Detailed dis- 
cussion of the Selb erg trace formula can be found in 
(lHeihall.ll976lll983h or mo re physics oriented approaches 
cmannl 



in ( Stockmaiml 19991) and (Wardlaw and Jaworski, 



1989) 



1993; 



and in the context of the Casimir - effect in ( Elizalde 
Kurokawa and WakavamaL 2002 : SchadenL l2006h . 

So, let us return to the scattering formalism in the 
standard Euc lidean space. Joffily , moti vated by Pavlov 
and Fadeev ( Pavlov and Fadeev . 19751) . examined the 
scattering states of a non-relativistic, spinless particle 
under the influence of a spherically symmetric, local 
and finite potential. He e xamine d the Jost solutions 
of this scattering problem ( JoffilvL 12003 ). which differ 
from the physical sol ution of the Schrodinger equation 
in th eir asymptotical (jAlfaro and Regge , 19651: iNewton , 
19821 ). In standard non-relativistic scattering theory the 
S-matrix is given by 



S(k) 



_ 2iS(k) _ 



f-(k) 
f+(k) 



(29) 



where S(k) is the phase shift, and f±(k) are the 
Jost solutions defined by t heir boundary co n dition s 
lim^oo (/±(fc)e =Flfcr ) = 1 (jAlfaro and Reggel . Il965l) . 
Provided the potential has a finite range and decreases 
sufficiently rapidly, the Jost solution f+(k) is proven to 
have infinitely many zeros, corresponding to the solu- 
tions of the Schrodinger equation as outgoing or incom- 
ing waves. Resonances (i.e. states with finite lifetime) 
occur if S(k) has poles on the complex k plane with neg- 
ative imaginary parts: k 2 — e n ~ iT n /2, where e„ and 
r„ stand for the energy and inverse lifetime associated 
with the nth state. Joffily introduces a mapping between 
these zeros of f+(k) onto the critical line and shows they 
coincide with the non-trivial zeros of the Riemann zeta 
function. He associates this artificial system with a vac- 
uum and the zeros are interpreted as an infinity of vir- 
tual resonan ces, and thus refle ct the chaotic nature of 
the vacuum (|joffilvl . [2003ll2004l) . This interpretat ion has 
also been extended using relativistic scattering (jjoffilvl . 
20071) . 

In another scattering based approach, Chadan and 
Musette analysed the so-called 'coupling constant spec- 



The Jost functions are the solutions of the Schrodinger equation 
with the following asymptotic behaviour: 



lim (e lkx f(X, k,x)) = 1 



where A : 



i is the shifted angular momentum, k ~ \[E and 
x £ (—00,00). This choice of the boundary condition is moti- 
vated by our physical picture, i.e. the particle should be repre- 
sented by free plane waves far from the local potential. The real 
physical solution of the Schrodinger equation can be expressed 
as a linear combination of the two Jost functions. 




Re(C(1/2 + it) 



FIG. 6 Argand diagram of the Riemann zeta function on the 
critical line, £(l/2 + ii), where t = 0—49.77, the later of which 
is £s pio. 



trum' o f a radially symmetric three -dimensional Hamil- 
tonian ( Chadan and Musette! . Il993 ) where the potential 
is chosen from a singular family of functions 



T~Lcm — — — 



dr 2 



CM 



(30) 



where /cm has logarithmic singularities at r = 0. They 
argued that the coupling constant spectrum coincides 
'approximately' with the non-trivial Riemann zeros if 
the problem is restricted to a finite, closed interval 
r € [0, e~ 47r / 3 ]. Mathematically rigorous detailed anal- 
ysis and the ex t ension of the potential family was carried 
out by (|Khuril . 2002 ). Furthermore, the existence of a 
three dimensional potential, Un(r), was derived whose s- 
wave scattering amplitude has the complex zeros of the 
Riemann zeta function as "redundant poles" . Exami- 
nation of C(s) using a quantum scattering approach is 
further motivated if one compares the plot of the phase 
of £(s) on the complex plane with with the usual Argand- 
diagrarr@ of the scattering amplitude corresponding to a 
collision. 

As a specific example, for completely elastic collisions, 
the scattering amplitude should be a perfect circle on the 
complex plane with unit radius centred on (0,1). For in- 
elastic collisions this circle deforms. The phase of ((s), 
after interchanging the roles of the real and imaginary 
axes, qualitatively resembles the Argand diagram of a 
scattering amplitude. This geometric similarity suggests 
an analysis of £(s) as if it represented the scattering am- 
plitude of a real collision of particles. This analogy, how- 
ever, is not perfect since £(s) does become negative while 



9 Argand diagrams can be thought of as a parametric plot of the 
inherently complex scattering amplitude on the complex plane, 
and the collisio n energy plays the role of the parame ter (see for 
example llBohm and Loewd . I2OO1T ) or llBhaduril. fl988t) ). 
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the Argand diagram of the scattering amplitude corre- 
sponding to a realistic collision does not. Bhaduri ad- 
vocates neglecting these small differences which do not 
affect their most important result, namely the phase 
9(t) of the Ricmann zeta function along the critical line, 
£(1/2 + it) = Z{t)e~ ie ^\ is intimately connected to the 
quan tum scattering of a particle on a saddle-like sur- 
face (IBhaduri et al. 



g ot a partic 
a ll995lll997t) 



To illustrate this, let us consider a non-relativistic par- 
ticle moving in an inverted harmonic oscillator potential 
along the half-line (x > 0). The Schrodinger equation 
reads as 



2m dx 2 



1 



$(a;) - -moj x $ (x) = E${x) 



(31) 



where we require that <$>(x = 0) = 0. This problem can 
be mapped onto a repulsive Coulomb problem of which 
the p hase shift 8{t) can be exactly expressed ( Fliiggel . 
19741 ) . The oscillatory part of the phase shift is given by 



8{t) = Smooth (*) + 



]n[r|i + 4 



In r 



(32) 



which is exactly the phase of the Riemann zeta func- 
tion. Two years after their first result, Bhaduri et al. ex- 
tended this one-dimensional model to a two-dimensional 
one where in one direction the potential is a traditional 
confining parabolic potential, and in the perpendicular 
dire ction (y) they kept t he inverted harmonic oscilla- 



tor (|Bhaduri et all 119971) . This choice was motivated 



by the analysis of the Gutzwillcr trace formula on the 
a = 1 border of the critical line, and also by the form of 
the electrostatic potential at the bo ttleneck of a qu antum 
contact in a mesoscopic structure ( Biittiker . 1990t ). 

While the inverted oscillator reproduced the oscillat- 
ing part of the £(s) phase in Bhaduri's work, Berry 
and Keating showed that a regularisation of a surpris- 
ingly simple one-dimensional classical Hamiltonian, % = 
xp, reproduces the smooth cou nting function of the ze- 
ros ( Berry and Keatingl Il999a ). We note here that this 
choice of T~L is a canonically rotated form of the inverted 
oscillator Hamiltonian ~ [p 2 — x 2 ). Moreover, the quan- 
tum mechanical model of the corresponding symmetrised 
Hamiltonian, % = (xp+px)/2 has also been investigated 
and exactly solved preserving the self-adjoint property 
of the Hamiltonian ( SierraL 120071 Twamlev and Milburnl . 
20061 ). The beauty of the xp- or inverted oscillator model 
is that it satisfies most of the properties listed earlier (see 
page I!]): valid as a classical mechanical model; the dy- 
namics is one-dimensional and uniformly unstable since 
the solution of the Hamiltonian equations are exponen- 
tially decaying or diverging; it lacks time-reversal symme- 
try. However, the trajectories are not bounded causing 
significant hardship in the semiclassical quantisation. As 
in the hyperbolic case, the boundary conditions or the 



way the phase space is regularised/ compactified become 
decisiv e. Berry and Keating sugge sted a simple regulari- 
sation ( Berry and Keatingl 1 1999al ) by introducing a cut- 
off in both position and momentum. This process results 
in a finite area, which can be filled up with Planck-cells 
of size h, thus counting the number of available quan- 
tum states. Another approach is available if one notices 
the dilation symmetry of the Hamiltonian (x i— > Ax, and 
p I—)- p/A). This symmetry manifests itself in the trans- 
formation of the wavefunction as 



ifj(Xx) 



Ai/2- 



Tlp(x) 



(33) 



and one might suggest restricting ourselves to A being 
a positive integer. This could be an attractive sugges- 
tion, because the wave-packet, generated by the uniform 
superpositions of all these transformed wavefunctions is 



¥(z) = X>(Az)=C(--i£U(s) 



A=l 



1 



(34) 



However, there is no physical motivation which would 
require this £ pre-factor to vanish. Furthermore this 
integer-based dilation-symmetry does not form a group, 
because the multiplicative inverse element (which would 
be A = 1 jm) is missing. 

Berry and Keating also established a peculiar canonical 
transformation (X = 2ir/p, P = xp 2 /2ir) for this Hamil- 
tonian, which exchanges and mixes the roles of the phys- 
ical position and momentum, but was uncertain "how to 
convert this "qu antum exchange" into an eff ective bound- 
ary condition" ( Berry and Keatingl Il999b| ). Aneva also 
analyses this boundary condition for a hyperbolic dynam- 
ical system with conformal geometry and shows how this 
exchange t ransformation arises as a result of boundary 
conditions (|Aneval . ll999L l2001al|bh . 

Later Sierra generalised Berry's model in two differ- 
ent ways, first by incorporating the fluctuation terms, 
d osc (E) (see re l ation (19b)) ). via changed boundary con- 
ditions (jSierral 120081 ). Together with Townsend, they 
also considered the motion of a charged particle (elec- 
tron) moving on the [xy] plane in a constant uniform 
perpendicular magnetic field, and in an electric potential 
described by the following Hamiltonian 



2/i 



pI 



eB 

Py + —X 



eXxy. 



(35) 



In this model, the number of semiclassical quantum 
states with energy less than E has the same functional 
form as the counting function of the Q(s) zeros, eq.[TUl i.e. 
the smooth part of the Riemann-zeros is reconstructed 
by the lowest-lying Landau level of the charged parti- 
cle. The fluctuation term - as they speculate - might be 
explained by the contribution of higher Landau-levels. 
This surmise, however, is only supported by estimat- 
ing the order of magnitude of these higher contribution 
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and comparing it to that of the Riemann £(s). This 
model has the additional attraction of being potentially 
accessible to ex perimentalists, including in lower spa- 
tial dimension s (|Li and Andrei , 2007 ; Park et. all l2009t 



Toet et. all 11991 ) 



Exploiting the x -H- p exchange symmetry of this model 
and using the Riemann-Siegel formula for the £(s) func- 
tion, Sierra created a new model in which the Jost solu- 
tions are directly proportional to the Riemann zeta func- 
tion, and the non-trivial zeros become the energies of the 
bound states. This achievement does not, however, prove 
the Riemann Hypothesis, as Sierra explicitly states "we 
cannot exclude the existence of zeros outside the critical 
line". 

In summary, we first introduced, motivated by 
Gutzwiller's trace formula, a quantum mechanical model 
on a surface with negative curvature, which lead us to the 
mathematically exact Selberg trace formula. The impor- 
tance of this result is at least twofold. Firstly, it reassures 
us that describing chaotic systems via the periodic orbits 
is likely to be feasible, and secondly demonstrates the 
role of periodic orbits in a generic system in determining 
the smooth and fluctuating parts of the density of states. 
We further elaborated on another non-Euclidean model, 
proposed by Pavlov and Fadeev, in which the Riemann- 
zeta function determines the 5-matrix over the complex 
energy plane. 

Converting these results into the usual Euclidean 
space, however, seems challenging. Although a few mod- 
els have successfully reproduced the smooth part of the 
density of quantum states, the fluctuation terms of these 
models differ from that of the Riemann zeta function. 



2. Bound state models 

From the 1950s a new approach, the Random Ma- 
trix Theory, emerged from the study of the spectrum of 
heavy nuclei. The same statistical apparatus had also 
been used to analyse the statistical properties of the 
seemingly random Riemann zeta zeros, and lead to the 
conjecture that th e ((s) zeros belong to one p articular 
universality class (jBohigas et al. , 1984a. 1986), the so- 
called Gaussian Unitary Ensemble (see later in section 
IHI.Cf) . This result suggested property [3] on Berry's list. 
However, Wu and Sprung generated a one-dimensional, 
therefore integrable, quantum mechanical model which 
can possess the Riemann ze ta zeros as energy eigenval- 
ues (|Wu and Sprung] . Il993l ) and show the same level- 
repulsion as that observed in quantum chaos. This was 
a contradictory result since on one hand the £(s) zeros 
follow a statistics specific for systems violating time re- 
versal symmetry, on the other hand, Wu and Sprung's 
model, by definition, was invariant under time reversal. 
However, the proposed model was not lacking in irregu- 
larity, since the potential reproducing the Riemann zeta 



zeros appeared to be a fractal, a self-similar mathemat- 
ical object. Nevertheless, these authors derived for the 
first time a smooth, semi-classical potential which gener- 
ate the smooth part of through 



N(E) 



dx dp = 



y/E - V(x) dx (36) 



with the 2m /h 2 set to unity. Solving this Abel- type in- 
tegral equation one may derive the following implicit ex- 
pression for V{x) 



X(V): 



x/V-V \n^ 



Fin 



V + y/V-Vo 



V-y/V-VoJ 

(37) ' 

where Vo has to be chosen such that the potential is not 
multi-valued, i.e. Vo < 2tt. The choice of Vo affects the 
potential at its bottom (x ~ 0), but for large x it does 
not have a significant impact and for i»l 



x(V) 



7r V 7re z 



(38) 



(see figure 1 in ( Wu and Sprung] . 1993 )1. We note here 
that Mussardo, using similar semiclassical arguments as 
Wu and Sprung, recently also gave a simple expres- 
sion f or a smoo t h pot ential supporting the prime num- 
bers (jMussardol . Il997f) as energy eigenvalues. Further- 
more, Mussardo also proposed a hypothetical resonance 
experiment to carry out primality testing; this theoreti- 
cally infinite potential could be truncated at some high 
energy; thus transforming it into a finite well. If an in- 
cident wave radiated onto this well has energy E = nfku 
where n is a prime number, then it should cause a sharp 
resonance peak in the transmission spectrum - argues 
Mussardo. 

Turning back to the smooth potential studied by Wu 
and Sprung, which is able to 'roughly' reproduce N(E), it 
is then modified to have the low lying £(s) zeros exactly. 
In order to achieve this goal Wu and Sprung set up a 
least-square minimisation routine, to minimise the differ- 
ence between the actual energy eigenvalues and the exact 
zeros. The result was surprising, since the potential curve 
became coarse and resembled a random potential. They 
analysed this curve using the standard box-counting tech- 
nique and measured a d = 1.5 fractal dimension for the 
potential reconstructing the Riem ann zeta zeros. 

Ramani et al. pointed out (jRamani et all [1995) 



that the apparent contradiction between Berry's conjec- 
ture and Wu and Sprung's model, i.e. whether or not 
the physical system exhibits time-reversal symmetry, is 
caused by the coarse curve of the potential, since any 



10 In the mathematical literature the argument is usually denoted 
by T, as in section [IT] Motivated by physics, we use here E. 
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smooth one-dimensional potential would lead to locally 
evenly spread energy levels, which is not the case for 
Wu's potential. They also provided a very efficient algo- 
rithm, the "dressing transformation" with which one can 
build up the quantum potential from individual energy 
eigenvalues. However, they standardised the spectrum 
using the 'spectrum unfolding' technique which eventu- 
ally lead them to the conclusion: the fractal dimension 
of the potential supporting the Riemann zeta zeros has 
d — > 2 r ather than that measur ed by Wu and Sprung. In 
a reply ( Wu and Sprune . 19951) . Wu and Sprung pointed 
out that this difference in fractal dimension is putatively 
caused by the alternative choice of spectrum. As they ar- 
gued, Ramani's spectrum does not have the same average 
density, long range correlation and nearest-level spacing 
distribution as the Riemann zeta function, therefore one 
cannot draw valuable conclusions regarding the potential. 

Nearly a decade after Wu and Sprung's origi- 
nal article, van Zyl and Hutchinson attempted to 
clarify the questions raised by th e two previous 



works ( van Zyl and Hutchinson . 20031) . They showed 



that for the same set of energy levels different poten- 
tial generating techniques (the variational approach used 
by Wu and Sprung, the dressing-transformation used by 
Ramani et al.) lead to the same potential, depicted in 
Figure [7] This result had been further strengthened by 
Schumayer et al. who used the inverse scattering trans- 
formation as a third techni que obtaining the s ame p oten- 



tials as in the earlier works (jSchumaver et all 12008). It is 



noteworthy to mention that the inverse scattering trans- 
form guarantees the uniqueness of the potential in one- 
dimension. This analysis, therefore, elucidated that the 
difference in measured fractal dimension cannot originate 
from the method of inversion. Moreover, they confirmed 
d = 1.5 for the Riemann zeta potential. These works all 
demonstrated the importance of long-range correlations 
in determining the fractal dimension of the potential. 

In a similar manner to that of the the Riemann £(s) 
zeros, the prime numbers can also be considered as an 
energy spectrum, thus a potential can be associated with 
them and it also proves to be fractal, but with a larger 
fractal dimension, d = 1.8. This result is somewhat puz- 
zling. The two sets, those of the zeta zeros and the prime 
numbers can be mapped onto each other via eq. (fB"]). 
but the nearest-neighbour spacing distribution of prime 
numbers is known to be Poisson-like (almost uncorrelated 
random distribution) while that of the Riemann zeros is 
rooted in the Gaussian Unitary Ensemble, and exhibits 
the corresponding correlations (see expression (|40| in sec- 
tion llll.Cl) . One may, therefore, conclude that Riemann's 
formulae converts two very different random distributions 
into each other, or as Sakhr et al. put it (jSakhr et all 
20031 ): "it is possible to generate the almost uncorre- 



lated sequence of the primes from the interference of the 
highly-correlated Riemann zeros" . 

Regarding the fractal nature, Schumayer et al. also 




nversion potential with N=200 - 
Semiclassical approximation ■ 
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FIG. 7 Main figure shows the semi-classical potential (dashed 
line), and the fractal potential (solid line) supporting the 
first two hundred zeros of ((s) as energy eigenvalues. The 
inset depicts the diffe rence of these potentials. From 
(jSchumaver et aUl200Sl ). 



established that the potentials associated with either 
the zeros of £(s) or with the prime numbers are multi- 
fractals, i.e. these potential curves cannot be charac- 
terised by one number d, but a range of dimension is 
necessary to describe th eir properties (for definition see 



(jSchumaver et all [2008)) 



Finally, at the end of this section devoted to the quan- 
tum mechanical models of the Riemann zeta function, we 
briefly refer to another alternative spectral interpr etation 
of the zeros proposed by Connes ( Connesl . 19991) . Dur- 
ing the comparison of the Gutzwiller's trace formula for 
quantum mechanical systems and that of the £(s) func- 
tion we noticed the overall sign difference in d osc (see the 
negative sign in equation (|9bl) in front of the summation) , 
i.e. the contribution of the periodic orbits should be sub- 
tracte d and not ad ded to the smooth density of states, 
d(T) dBerrvl . Il986h . This sign difference led Connes to 



interpret the zeros as gaps, missing lines from the oth- 
erwise continuous energy spectrum rather than discrete 
energy levels. 



C. Nuclear physics 

Random Matrix Theory (RMT) has been successfully 
applied to predict ensemble averages of observables for 
heavy nuclei. Even though the Riemann zeros are 
distributed randomly, some of their statistical quanti- 
ties correspond to that of the Gauss Unitary Ensem- 
ble. We discuss the RMT briefly for historical rea- 
sons. The reason for brevity owes to two recent Co l loquia 
devoted to RMT IPavenbrock and Weidenmuller . 200$ : 
Weidenmuller and Mitchell \200A) . 

Unfortunately the degrees of freedom of even a mod- 
erately large nucleus are still far beyond our computa- 
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tional capability, be it analytical or numerical. Similar 
problems, although the number of components are on a 
different scale, have occurred before in physics and en- 
gendered the development of a new branch of physics, 
statistical mechanics. This is exactly what Wigner had 
in m ind when he su ggested a statistical description of nu- 
clei (|Wignerl . Il951l) . He suggested that nuclei can be sta- 
tistically described using random matrices carefully cho- 
sen from pre-determined ensembles. The new description 
emerging from this examination is the Random Matrix 
Theory. 

Although random matrix theory emerged from the sta- 
tistical description of nuclei, it has already infiltrated into 
many different areas of physics. Recent develo pments of 
this branch of physics have been reviewed in (Bohigas, 
19891: iForrester et all [200l IWeidenmuller and Mitchell. 



20091 ). Moreover we can suggest the monograph by one 



of th e leading figures of random matrix theory ([Mehta 



2004) 



But how to choose the ensemble of random matri- 
ces suitable for a certain system, or for the Riemann 
£(s) function? Throughout classical mechanics symme- 
try plays a decisive role in determining the dynamics of 
different systems. If a physical system has a symme- 
try it implies, via Noether's theorem, the existence of a 
conserved quantity, e.g. the translational invariance in 
time dictates energy conservation, continuous rotational 
invariance requires the angular momentum remain con- 
stant. These symmetries limit the possible forms of the 
Hamiltonian describing the given system. Therefore, if 
one wants to approximate a Hamiltonian with a large, 
but finite dimensional matrix these symmetries will de- 
termine the type and structure of the matrix, whether 
it is real or complex, symmetric or hermitian (jDvson , 



19G2) 



In the case of an integrable system, the conserved 
quantities are all known. Therefore the Hamiltonian can 
be diagonalised, with each eigenvalue forming its own 
symmetry-class. This leads to the assumption that these 
eigenvalues are completely uncorrelated. Let us also as- 
sume that the average spacing between eigenvalues is 
unity in the overall sequence of eigenvalues. If p(s) de- 
notes the probability distribution of nearest neighbour 
spacings, i.e. if ei and £2 are eigenvalues of t he given sys- 
tem, then ei — £2 = s, then one can express (jStockmann 



Il999h the probability of finding two eigenvalues in a dis- 
tance between s and s + ds with no other eigenvalues in 
between. Dividing the distance s into N equal intervals, 
the probability is simply 



p(s)ds — lim 



iV 



(39) 



In the N — > 00 limit the right hand side becomes the 
exponential function. Therefore the probability distri- 
bution p(s) = exp(— s) is the Poisson distribution with 
parameter equal to 1. This is quite a general result 



for integ rable systems as Berry a nd Tabor have demon- 
strated ( Berry and Tabori . Il977h . Similarly, one can 
deduce similar probability distributions for universality 
classes of random matrices, e.g. Gaussian Unitary En- 
semble, Gaussian Orthogonal Ensemble, etc. The classi- 
fication refers to the universality conjecture: if the clas- 
sical dynamics is integrable then p(s) corresponds to the 
Poisson ensemble, while in the chaotic case p(s) coincides 
with the corresponding quantity for the eigenvalues of 



a suitable ensemble of random matrices (Bohi gas et 



1984al) . Furthermore, the local statistics of the eigenval- 
ues converge as the order of the matrix increases. 

How is this connected to the Riemann-zeta zeros? 
The zeros can be treated as eigenvalues of a fictitious 
physical system, just as Hilbert and Polya suggested, 
and their statistical p roperties examined. In 1973 Hugh 
Montgomery showed (j Montgomery . Il973h that the pair- 
correlation function of the zeros is 



r 2 {x) = 1 - 



sin (ttx) 



(40) 



provided the Riemann Hypothesis is true. Free man J . 
Dyso n, during an informal discussion over tea ( Cipral . 
1999), pointed out to Montgomery that this is exactly 



the same result as obtained for random matrices picked 
from the Gaussian Unitary Ensemble. However, this 
statement is made in the asymptotic limit, i.e. as one 
goes to infinity on the critical line, 1/2 + iE, or in the 
RMT language as the size of the matrices, N, tends to 
infinity. At finite height E or dimensionality N discrep- 
ancies may occur compared to expression (|40[) . Inter- 
estingly, it was shown using heuristic arguments, that 
the nearest-neighbour spacing distribution of the zeta 
zeros and that of unita ry random matrices of finite di - 
mension are the same ( Bogomolnv and Keating . 19951) . 
Moreover, the sam e authors extended their study of cor- 
relation functions ( Bogomolnv and Keating . 19961 ). r n of 
order n (n > 2) and proved, in the appropriate asymp- 
totic limit, r n of the Riemann zeta zeros are equivalent 
to the corresponding GUE result. This result was com- 



plimentary to Montgomery's second order ([Montgomery , 
Il973h . Hejhal's third order (jHeihall . Il994h and Rudnick 
and Sarn ak's general result for the n th order correlation 
function (jRudnick and Sarnakl . Il996). 

What does this result demand from a model of the Rie- 
mann zeros? The striking similarity between the pair- 
correlation function of the £(s) zeros and the eigenval- 
ues of random matrices from the GUE ensemble only 
holds for short-range statistics. Odlyzko, by calculat- 
ing the statistics for subs tantial number s of zeros, car- 
ried out an empirical test riOdlvzkd . [l987h and confirmed 
Berrvl 



Berry's predictions (jBerr 

v. 1 19851) about the discrepan- 



cies between the GUE theory and computed behaviour 
of the £(s) zeros. The long-range correlation and the 
small spacing statistics of the ((s) zeros noticeabl y devi- 
ate from the GUE prediction. This is expected (jBerrv . 
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1985, 1988), since long-range correlations are dominated 
by the short periodic orbits, which are system specific 
and therefore not universal. For £(s) the mean separa- 
tion between zeros is ln(£'/27r) while the smallest period 
is ~ ln(2) (Table [TJ) . Conclusively, the GUE-predicted 
universal correla tion for zeros near E shoul d fail beyond 
ln(£/27r)/ln(2) (|Berrv and Keatinsi[l999b]) . Despite the 
deviation explained above, the statistics of the £(s) zeros 
asymptotically coincide with those of the GUE ensemble, 
consequently the corresponding qua ntum system ought 
to vio late time-reversal symmetry (jBerrv and Keating , 
1999al lbh. This may have motivated Berry and Keating's 
choice of a ~ (xp + px) as a Hamiltonian. 

Finally, we must mention an unexpected spin-off re- 
sult of random matrix theory related to the P olya-Hilbert 
conjecture. Crehan asserted ( Crehanl . Il995l ) that for any 
bounded sequence there are infinitely many classically in- 
tegrable Hamiltonians for which the corresponding quan- 
tum spectrum coincides with this sequence. Further- 
more, as an example for his theorem, he shows that in- 
finitely many classically integrable non-linear oscillators 
are capable of exactly reproducing the Riemann-zeta ze- 
ros when they are quantised. Unfortunately, the theo- 
rem is an existence theorem and not a constructive one. 
If such a system could be created, whether physically 
or just theoretically, that would be aesthetically pleas- 
ing: it would connect the most studied physical model 
(oscillator) with the basis of our arithmetic (prime num- 
bers). Crehan's result is promising and is also supported 
by the relationship between the Riemann £(s) zeros and 
the Painleve V equation, the latter of which plays a cen- 
tral role i n the theory of completely inte grable dynamical 
systems ( Ablowitz and Clarksonl . [l99lh . 

Finally, in this section we briefly mention the notion 
of quantum ergodicity which attracted substantial atten- 
tion in the last three decades in the search for links be- 
tween classical and quantum ergodicity, i.e. what "fin- 
gerprint" the classical chaos leaves in the physical prop- 
erties if we quantise the system, especiall y in the long- 
time behaviour. Only few rigorous results (Schnir ehnan , 
1974 Ide VerdiCTeUl985UZelditch and Zworskil . Il996l) are 
known, and one of them says that the expectation value 
of operators over individual cigcnstatcs is almost always 
the ergodic, microcanonical average of the classical ver- 
sion of the operator. However, the theoretically rigor- 
ous understanding of quantum ergodicity is still in its in- 
fancy. Numerical simulations suggest though that quan- 
tum chaotic systems exhibit universal behaviour at a 
particular length scale, and at this scale the statistics 
of the eigenvalues resemble that of l arge random matri- 



ces chosen f r om specific en s embles (lAgam et al. 



Berry 



Heller 



1977E iBohigas et ail ll984allbinGutzwiller 



1995 



1991 



1984). It is unfortunate that this length scale 



19911 ) that quantum eigenstates must deviate from the 
RMT predictions. These corrections may stand out from 
the spread out background of RMT, just as the unstable 
periodic orbits do as eigenstates with enhanced ampli- 
tudes as depicte d in Figure [51 Although further numer- 
ical s imulations ( Backer et all 11998c iKaplan and Heller , 
1999h provide some evidence regarding the connection 
between RMT and quantum ergodicity, its interpretation 
and strength remain open questions. 



D. Condensed matter physics 

In condensed matter physics the fundamental structure is 
the crystal lattice. Below we examine the connection of 
the lattice with the generalised Riemann hypothesis. We 
also show how the specific heat capacity of a solid restricts 
the location of the £(s) zeros. 

One of the fundamental bases of modern condensed 
matter physics is the geometrical structure of solids; the 
lattice. The examination of this mathematical struc- 
ture is necessary to understand even the basic proper- 
ties of matter. The regular structure of a perfect lat- 
tice is suitable for immediate comparison with regulari- 
ties among the natural numbers, and therefore it is not 
a surprise that many number-theoretical functions arise 
in crystallography, e.g. Ninham et al. present a witt y 



review on the Mobius function ( Ninham et al . 19921 ). 



is so minute that it hinders the numerical simulations 
substanti ally. Nevertheless, it has also been shown theo- 
retically ( Kaplan and Heller . 19961 : Tomsovic and Heller , 



For those mathematically more inclined we suggest the 
book "From Number Theo ry to Physics" by Waldschmidt 
(|Waldschmidt et aLl . [l995h . Moreover, not only the per- 
fect regularity of a lattice, but also the lack of this reg- 
ularity can be related to th e Riemann z eta function, 
as Dyson indicated recently ( Dvson . 2009f ): "A fourth 
joke of nature is a similarity in behavior between quasi- 
crystals and the zeros of the Riemann Zeta function." In 
the following, we briefly examine why a solid state con- 
stituted by ions should even exist, what binds these ions 
to each other? 

Ions arrange themselves into a structure which max- 
imises the attractive interaction between unlike and min- 
imises the repulsive interaction between like charges. In 
an ionic crystal, such as NaCl, the main contribution to 
the binding energy has an electrostatic origin with the 
van der Waals term only a few percent of the former. 
The electrostatic term is called the Madelung energy, and 
the energy of one ion in the solid is called the Madelung 
constant. 

For the sake of simplicity, let us first imagine a one- 
dimensional infinitely long ionic lattice. Cations and 
anions are located next to each other at a distance 
a, in a simplified NaCl structure. If simply two unit 
charges q were positioned at the same distance a, the 
electric potential energy of one of the charges would be 
U = q 2 /47reoo. In a solid each ion is in the field of all 
the remaining charges, both positive and negative. The 
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FIG. 8 Schematic structure of a fictitious one-dimensional 
solid built up by cations and anions, positioned in alternating 
pattern. 



as Cep(1/2, <5 mi „) where £ep is the Epstein zeta function 
(see below), and m,n =1, 2, 3. The second argument, 
<5 mj „, is determined by the type of the lattice, and in crys- 
tallography it is a quadratic form V given by the Gram 
matrix p mn = e m e„, where e m is the mth lattice vector. 
Therefore, for example, the Madelung constant for the 
body-centered cubic structure can be formally written as 



total electrostatic potential energy of one ion at position 
i in the lattice is therefore 



Ui = 



x ^ i (-i)i'-y 

Z-f. 47re \i-j\a 



1 q 2 



(-1)* 



47ren a ^—f k 



(41) 



where j runs over all lattice sites except i in the first sum- 
mation, and in the second we have changed the running 
variable to k = \i — j\. In a finite lattice we have 2N ions, 
but in (|4"Tj) each term belongs to two ions, therefore the 
total electrostatic potential energy of the finite lattice is 



Wtotal = I 2NUi = X- 



1 



Q 2 ^ (-1? 



Airen a 



k=£0 



(42) 



This form of Wtotai can be divided into three terms: X 
which guarantees the extensive nature of the energy, an 
energy factor, q 2 /Altera, and also a numerical factor de- 
pending only on the lattice structure. One sees directly 
that the influence of the lattice on the total electrostatic 
energy is comprised of an infinite sum. Since this energy 
term has to be negative in order to describe binding, we 
incorporate this sign into the Madelung constant chid as 



ckid 



= 2£ 



(-1) 



fe+i 



(43) 



k=l 



where the factor 2 appears because of the mirror- 
symmetry around the ith ion. The total energy can be 
written as U% tal — — oiyoXq 2 /Altera, which is negative if 
aiD > 0. 

Generalising the NaCl structure we examined above 
for the realistic three-dimensional case, one can write the 
Madelung constant for this lattice as 



= 

(i,k, 0^(0,0,0) 



(-1) 



i+j+fc+l 



d 2 +J 2 



jfc a )Va- 



(44) 



Although it is tempting to evaluate this summation by 
approximating the terms on concentric spheres centred at 
the reference ion (i = j = k = 0) and utilising the sym- 
metry, the resulting series, 6-12/V2+8/V3-.. . is diver- 
gent which is physically unsatisfactory. The convergence 
proper ties of such sums have been extensively investi- 
gated dBorwein et all 11985k IChaba and Pathria . 19751 



1771 ). The sum (j44|) is an alternating and con- 
ditionally convergent sum. The denominator of the sum- 
mand is a quadratic form, therefore the Madelung con- 
stant for a simple cubic structure can be formally written 



n bcc 



-,EP 




= 1.762675. 



(45) 



Here we have only dealt with the pure Coulomb- 
interaction, but this treatment can be ex- 
tended to screened electrostatic interactions as 



well (jKanemitsu et ali , 120051 ) 



The infinite sum in (|4"3"1) strongly resembles the 
Riemann-zeta function, except each term is weighted 
by a factor (— l) k+1 , and its numerical value is a-yn = 
2 In (2) w 1.3863. Although in two and three dimensions 
the summation can be written explicitly, obtaining a pre- 
cise numerical value is far from easy and the Epstein-zeta 
function is required. This fu nction can be thought of as 
a generalised zeta- function ( Ivic . 20031 : Shanker, 2006 ) 
which is defined by 



Cep(s,P) 



E- 



(46) 



where V is a quadratic form defined on a d-dimensional 
lattice. All lattice points for which V = are excluded 
from the summation. This function can be analytically 
continued to the same domain as the Riemann-zeta func- 
tion and also has its only pole at s = 1 with residue 
7r/vA- The similarity goes further since Cep(s,V) also 
satisfies a functional equation expressing mirror symme- 
try. Thus, there is an inclination to generalise Riemann's 
conjecture: all non-trivial zeros of Cep(s,7-') have real 
part one half. The temptation to do so is strengthened if 
one chooses specific quadratic forms, e.g. Cep(s, Idi) = 
2£(2s), or CEp(s,Id/t) ~ ((s)((s — 1), where Id„ is the n- 
dimensional identity matrix. Indeed, it was shown eighty 
years ago that for binary quadratic forms (two dimen- 
sional lattic e), infinitely many zero s of Ce p lie on the 
critical line ( Potter and Titchmarsh . 1935f) in a similar 
mann er to Hardy for the Riemann zeta function ((Hardy, 
19141 ). Remarkably, however, it has also been shown 
that in any dimension one can construct such a V , tha t 
the generalised hypothesis does not hold ( Terras . 198fjf) . 
This, admittedly negative, result shows the intriguing 
connection between crystallography and this generalised 
Riemann Hypothesis, but let us now depart from the ab- 
stract and static crystal structure of solids, and examine 
the dynamics of this system. 

The lattice vibrations, phonons, are bosonic quasipar- 
ticles. Therefore if one knows their energy spectrum, 
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hujk, then the total energy of the phonon gas is simply 
the sum over all modes of the crystal 



C/ = E 



e hu} k /k B T _ i ' 



(47) 



Since the number of possible modes is large, 3iV, where 
N is the number of atoms in the lattice, one might con- 
vert this expression into an integral by introducing the 
phonon density of states g{io) normalised as J g(uj)dui = 
3N. Using standard methods to calculate the specific 
heat of the solid, a directly measurable quantity, the fol- 
lowing expression can be obtained 



cv 



hw/k B T 



k B T J (e^/ k BT _ i); 



g(uj)duj. (48) 



The only sample-specific quantity here is g(u>). Surpris- 
ingly, the number theoretical Mobius function and the 
related Mobius inversion provides a transformation to ex- 
press g(j-L>) as a function of the measured specific heat 



7H 



1 



cy(h/k B u) 



(49) 



where the inverse Laplace-transform, £ , converts the 



space of u = h/k B T to w/n (IChenl . Il990h Around 
the same time another inversion technique appeared in 
the literature ( Xianxi et all Il99fjf) an d was proven to b e 



equivalent to the one discussed above (jMing et all 120031) 



In the latter formulation another special function, the 
Riemann zeta function was used, but in order to avoid 
the dependence on the unproven Riemann Hypothesis a 
free "regularisation" parameter s was also introduced. 
The density of states in this formalism is 



Zttuj 



where Q(k) = 



.ik+s 



Q(k) 



T(ik 



2) C(ik 



1) 



dk (50) 



.ik+s — 1 



cv(l/u) du. Physically the 
density of states, g{uj) should be independent of the reg- 
ularisation parameter, although the existence of Q(k) re- 
quires that s must fall into the < s± < s < S2 range 
where Sj and S2 are the exponents of the specific heat 
asymptotes at high- and low temperatures, respectively. 
Due to the Dulong-Petit law, at high temperature the 
specific heat is independent of the temperature, therefore 
si = 0. On the other end of the temperature scale the 
specific heat of phonons vanishes as T d in d dimensions. 
Therefore C(s) in the denominator of the integrand ([50)1 
sweeps through the [1, 1 + d] strip and ensures that no ze- 
ros of C(s) can occur there. Summarising, the asymptotes 
of the specific heat contribution of lattice vibrations in a 
solid provide an experimentally determined zero-free re- 
gion of C(s) on the complex s plane. Although this offers 
no further restriction than that which is already known 
from mathematics, it is an example where physics places 
independent bounds upon the location of the zeros. 



E. Statistical physics 

The description of both bosons and fermions relies on the 
mathematical properties of the Riemann zeta function. 
We show how the problem of the 'grand canonical catas- 
trophe ' of number fluctuation in an ideal Bose-Einstein 
condensate is connected to number theory. We introduce 
the concept of the primon gas, and also consider number 
theoretical models of Brownian motion. 

Although statistical physics, the physics of systems 
with a large number of degrees of freedom, relied heav- 
ily upon combinatorics well before the birth of quantum 
mechanics, probably the first appearance of the Riemann 
zeta function in statistical physics occurred in Planck's 
momentous work on black body radiation, the dawn of 
the quantum era. From then on, the Riemann zeta func- 
tion pops up in numerous different branches of statistical 
physics, from Brownian motion to lattice gas models. 

Since the topic of ultra-cold quantum gases has ex- 
panded rapidly in the past decade, we interpret the im- 
plications of the distribution of the Riemann zeta zeros 
in this area first. We start with the non-relativistic, non- 
interacting, spin zero Bose gas and treat the spatial di- 
mension D as a free para meter. It is a standard textbook 
derivation ( Huand . 2001 ) to show that this system under- 
goes a phase transition at low temperatures, where the 
de-Broglie wavelength A of the particles becomes compa- 
rable to the inter-particle distance, and thus the quantum 
nature of the constituents becomes decisive. Since the 
particles are free, their spectrum is continuous, simply 
equal to the kinetic energy e = p 2 /2m. The total num- 
ber of particles N is the sum of particles in each quantum 
state 
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(51) 



Changing the integration from momentum to energy 
leads directly to 



N cx 



c D/2-l 



,(e-n)/kT 



— de oc CI — 
1 H 2 



(52) 



In the last step we used the fact that the chemical poten- 
tial approaches the energy of the lowest lying state, i.e. 
H = 0. 

This result shows that the Bose-Einstein condensa- 
tion phase transition cannot occur in homogeneous non- 
interacting systems in dimensions lower than three. The 
total number of atoms is a positive number and fixed for 
our system. In one spatial dimension, since £(1/2) < 
0, the positivity of N cannot be fulfilled. For two- 
dimensions the right hand side of ([52"T) is divergent due 
to the pole of the Riemann zeta function C(s) at s = 1, 
therefore N appears to be infinite. The position of 
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this pole can be interpreted as the manifestation of 
the Mermin-Wagner-Hohenberg theorem, which guaran- 
tees that a homogeneous two-dimensional system, pro- 
vided the interaction is sufficiently weak, cannot un- 
dergo a phase transition. One may thus see that the 
pole structure of the Riemann zeta function determines 
whether our system of interest can undergo a phase tran- 
sition or not. We note here that this phase transition 
can o ccur in lower dimensions for i nhomogeneous sys- 



Dai and Xid . 120031: 



terns (Baana to and Kleppnerl . 11991 
IWidoml . ll968l ). 

Let us turn to another fundamental question of sta- 
tistical mechanics: the equivalency of different statistical 
ensembles. The difference between the predictions for 
the "Riemann gas" (see below) based on microcanonical, 
canonical, and grand-canonica l ensembles has been inves - 
tigated by Tran and Bhaduri (jTran and Bhaduril 120031) . 
However, the motivation for the analysis is rooted in the 
so-c alled "grand-can onical catastrophe" of an ideal Bosc 



gas (jZiff et aLl . ll977t ). The number fluctuation of an ideal 



boson gas is 



(53) 



k=0 



where (nk) denotes the ensemble average of the occupa- 
tion number of the fcth energy eigenstate. According to 
the formula above, in the presence of a macroscopically 
occupied ground state, the number fluctuation is propor- 
tional to the total number of particles, SNq ~ N, which, 
in the thermodynamical limit (N — > oo), leads to diver- 
gence. 

Grossmann and Holthaus examined the illustra- 
tive model system of an ideal Bose gas trapped 
in a d-dimensional potential with a power-law en- 



ergy spectrum, e L 
energy eigenstates 



hbJl 



where 



lab els the 

( Grossmann and Holthausl . 1997b; 
I Weiss and Wilkens . 19971 ) . Later Eckhardt extended the 
analysis to the mean density of states and the level 
spaci ng distribution for ideal quantum gases (jEckhardt , 
1999). Grossmann et al. showed how the dimensional- 
ity and a, which, in some sense, measures the strength 
of the potential, depress or enhance the number fluctu- 
ation of the ground state as a function of the rescaled 
temperature t = ksT '/frio: 



(SN f 



Ct d l° 
t 2 ln(i) 
C{d/a-\)t 2 



(0 < d/a < 2) 
(d/a = 2) 
(2 < d/a) 



where C is calculated from a d- dimensional Epstein zeta 
function ([Holthaus et all 120011 ) , although here its value 
does not play a significant role. Therefore in a given 
spatial dimension the potential can enhance the fluctu- 
ation while dimensionality depresses it. They also ex- 
amined the behaviour of the heat capacity around the 



critical temperature to and proved that the heat capac- 
ity changes continuously at to if 1 < d/a < 2, but if 
d/a > 2 it undergoes a jump given by 



Nk f 



con 



c(f-i) 



(54) 



where C< and C> denote the asymptotic values of the 
heat capacity at t — > to from below and above, respec- 
tively. It is worthwhile to note that (5 No) 2 in the canon- 
ical ensemble could be expressed as the following integral 
over the complex plane 
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(t)A(/M)C(* - 1) 



(55) 



where A((3,t) = ^ (/3e n )~* is the spectral zeta function 
of a given spectrum e„ and r is chosen so all the poles 
of the integrand lie on the left of the path of integra- 
tion. Therefore all the results shown above are deter- 
mined by the pole structure of the spectral and the Rie- 
mann zeta functions, A(/3,£) and C,(s), respectively. The 
large-system behavior is extracted from the leading pole, 
while the finite-size corrections are encoded in the next- 
to-leading poles. 

The formulae (|54[) and ([55]) above did not just 
clarify an important physical question, namely num- 
ber fluctuation properties of a d-dimensional bo- 
son gas below the critical temperature, but also 
had valuable number theoretical consequences. The 
probl em solved above is a purel y combinatoria l 
one ( Grossmann and Holthausl 1997al Holthaus et al. 



2001 



Weiss et all 120031 : IWeiss and Holthausl l2002h : how 



many ways can one distribute n excitation quanta over 
N particles? This question, for general n and N, is quite 
difficult. However, in the low temperature limit the num- 
ber of excitations, n, is much smaller than the number of 
particles, N. This problem thus becomes tractable and 
one could obtain the results mentioned above. Calcu- 
lating the number fluctuation of a boson gas in a one- 
dimensional (d = 1) harmonic potential (a — 1) pro- 
vides (SNq) 2 ~ t. But t is simply proportional to the 
number of energy quanta 'stored' in the excited states, 



t = (k B T/hu) = n and therefore (SN y 



A mathe- 



matician - according to Grossmann and Holthaus - can 
now interpret this formula: 

If one considers all unrestricted partitions of 
the integer n into positive, integer summands, 
and asks for the root-mean-square fluctuation 
of the number of summands, then the answer 
is (asymptotically) just yjn. 

An intriguing consequence of this analysis is that a Bose- 
Einstein condensate could be used (in theory at least) 
to factorise numbers ( Weiss et a/.ll2004h which could be 
treated as a quantum computer calculating the prime 
factors. 
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Furthermore, using their physical insight, Weiss and 
collaborators could derive the following non-trivial num- 
ber theoretical result. Let <fr(n, M) denote the number 
of partitions of n into M summands regardless of their 
order (e.g. $(5,2) = 2 while $(5,4) = 1), and Q(n) 
stand for the total number of different partitions, i.e. 
Q{ n ) — ELi TO )- ^ * s a natural step to intro- 
duce the "probability" of having exactly M terms in a 
random partition by p(n,M) — $(n, M)/fi(n). It was 
then shown that this probability distribution does not 
become Gaussian, and it adopts its limiting distribution 
shape if n > 10 10 , which itself is a remarkable fact. 

Here we only mention that the same combinatorial 
problem arises in many different branches of math- 
ematical physi cs, such as lattice animals in statis- 



2001: Wu et al 



tical physics (jLima and de Menezes 
1996), numerical analysis on combinatorial o ptimi 



sation (Andreas and Beichl. 120031: Bauke et al . 2003L 



Maiumdar and Krapivskvll2002t Mertensl . 



19981 ) and also 

in the description of the low-energy excitations of a one 
dimensional fermion -system as bosonic degree s of fr ee- 
dom (bosonisation) ( Schonhammer and Meden . 1996h . 

Tran and Bhaduri's, and then Holthaus and Weiss' 
works further underline that the irregular behaviour of 
the canonical ensemble lies in the combinatorics of parti- 
tioning integers and the microcanonical and canonical en- 
sembles prognosticate dramatically different ground state 
number-fluctuations, Suq. This is an important example 
which unequivocally shows that the standard statistical 
ensembles can not always be regarded as equivalent. 

These examples, while not directly related to any at- 
tempt to prove the Riemann Hypothesis, but rather just 
the zeta function, do illustrate that results in physics can 
have profound implications for mathematics in general 
and number theory in particular. 

The interpretation of prime numbers or the Riemann 
zeta zeros as energy eigenvalues of particles appears not 
just in quantum mechanics but also in statistical me- 
chanics. Below, we review two concepts: the Riemann 
gas, sometimes called the primon gas, and the Riemann 
liquid, although their definitions vary slightly. 

In 1990 Julia propo sed the idea of a fictitious, non- 
interacting boson gas (|.Tulial . Il990l ). where a single par- 
ticle may have discrete energy equal to eo, ei,. . .where 
e n = eoln(p„) (n > 1) and p n stands for the nth prime 
number. This is why the constituents are called primons. 
Since the particles are not interacting, a many-body 
state, in the second quantised formalism, can be repre- 
sented by an integer number n. This natural number has 
a unique factorisation, n — p™ 1 p™ 2 ■ ■ ■p k nk which tells 
us that mi particles are in the \p\) state, ni2 particles 
are in the \p2) state and so on. Due to this uniqueness, 
each many-body state is enumerated once and only once. 
Therefore, the total energy of the system, in the state \n), 

is E n = TOie ln(pi) + m 2 e ln(p 2 ) H h m k e ln(p k ) = 

eo ln^™ 1 ^™ 2 • • •p™ fc ) = eo ln(n). In order to describe this 



gas we have to construct the partition function from this 
spectrum 



^exp 

n=l 



E n 



oo 1 



(56) 



where s = eo/kBT — (3eo and (3 — (ksT^ 1 is the inverse 
temperature. The partition function for the primon gas 
is thus the Riemann zeta function £(s) and hence the al- 
ternative nomenclature. It is apparent, by looking at the 
domain of C(s), that Zb is well-behaving for s > 1, i.e. at 
low-temperatures, while s < 1 is physically unacceptable. 
The boundary, s = 1 represents a critical tempe rature, 
called the Hagedorn temperature ( Hagedornl . ll965l) above 
which the system cannot be heated up, since its energy 
becomes infinite 
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(57) 



A similar treatment can be built up for fermions rather 
than bosons, but here the Pauli exclusion principle has 
to be taken into account, i.e. two primons cannot occupy 
the same single particle state. Therefore rrii can be or 
1 for all i. As a consequence, the many-body states are 
labeled not by the natural numbers, but by the square- 
free numbers. These numbers are sieved from the natural 
numbers by the Mobius function. The calculation is a 
bit more complex, but the partition function for a non- 
interacting fermion primon gas reduces to the relatively 
simple form 



Z F = 



Of) 

C(2s) 



(58) 



The canonical ensemble is of course not the only en- 
semble used in statistical physics. Julia extends the 
study to the grand ca nonica l ense mble by introducing 
a chemical potential \x ( Julial . [l994 ) . therefore replacing 
the primes p with new 'primes' pe~^. This generalisa- 
tion of the Riemann gas is called the Beurling gas, af- 
ter the Swedish mathematician Beurling who generalised 
the notion of prime numbers. Examining a boson primon 
gas with fugacity —1 shows that its partition function is 
2£ = C(2*)/C(s). 

This last result has an astonishing interpretation. We 
know that for a system, formed by two sub-systems not 
interacting with each other, the overall partition function 
is simply the product of the individual partition func- 
tions of the subsystems. Equation (|58|) has precisely this 
structure, there are two decoupled systems: a fermionic 
"ghost" Riemann gas at zero chemical potential and a 
boson Riemann gas with energy-levels E n = 2eo ln(p n ). 

Julia also calculates the appropriate Hagedorn tem- 
peratures and analyses how the partition functions of 
two different number theoretical gases, the Riemann gas 
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and t he 'log gas' behave around the Hagedorn tempera- 
ture (jjulial . ll99J) . Although the divergence of the parti- 
tion function signals the breakdown of the canonical en- 
semble, Julia also claims that the continuation across or 
around this critical temperature can hel p understand cer - 



tain phase transitions in string theory ( Deo et al. . 1989) 



or in the study of quark confinement (|Julial . 11994). The 
Ricmann g mathematically tractable model, has 

been followed with much attention because the asymp- 
totic density of states grows exponentially, d(E) ~ e B , 
just as in string theory. Moreover, using arithmetic 
functions it is not extremely hard to define a transition 
between bosons and fermions by introducing an extra 
parameter, n which defines an imaginary particle, the 
non- interacting parafermions of order k. This extra pa- 
rameter counts how many parafermions can occupy the 
same state, i.e. the occupation number of any state falls 
into the [0, k — 1] range, thus k = 2 belongs to nor- 
mal fermions, while k — > oo represents normal bosons. 
The partition function of a free, non-interacting n- 
paraf ermion gas can be shown to be ( Bakas and Bowickl 
19911) 



Z K (s) 



CM' 



(59) 



Bakas further demonstrates, using the Dirichlet con- 
volution (★), how one can introduce free mixing of 
parafermions with different orders which do not interact 
with each other 



(60) 



d \ n 



where the shorthand notation d\n means d is a divisor of 
n. This operation preserves the multiplicative property 
of the classically defined partition functions: Z Kl + K2 = 
2 Kl Z K2 . It is even more intriguing how interaction can be 
incorporated into the mixing by modifying the Dirichlet 
convolution with a kernel function or twisting factor 



.g = Y J f{d)g(- d )K{n,d). 

d\n 



(61) 



Using the unitary convolution Bakas establishes a ped- 
agogically illuminating case, the mixing of two identical 
boson Riemann gases. He shows that 
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(62) 



Thus mixing two identical boson Riemann gases inter- 
acting with each other through the unitary twisting, is 
equivalent to mixing a fermion Riemann gas with a bo- 
son Riemann gas which do not interact with each other. 
This leads to the interpretation that one of the original 
boson components suffers a transmutation into a fermion 



gas. It is noteworthy to mention that the Mobius func- 
tion, which is the identity function with respect to the * 
operation (i.e. free mixing) reappears in supersymmetric 
quantum field theories as a possible representation of the 
(— 1 ) F operator, where F is th e fermion number opera- 
tor jSpectorl . Ii~989l Il990l fl~998h . In this context, the fact 
that n(n) = for square-free numbers is the manifesta- 
tion of the Pauli exclusion principle. 

It is therefore interesting that, what initiated as rather 
academic studies to investigate potential attacks on the 
Riemann Hypothesis, may lead to advances in physics. 
But let us return to the Hypothesis through a slightly 
different definition of the Riemann gas. Here the energy 
of the ground state is taken to be zero and the energy 
spectrum of the excited state is e„ = ln(p„), where p n 
(n =2, 3, 5, . . . ) runs over the prime numbers. Let 
N and E denote the number of particles in the ground 
state and the total energy of the system, respectively. 
As we demonstrated above, the fundamental theorem of 
arithmetic allows only one excited state configuration for 
a given E = ln(rt) (n is an integer). It immediately 
means that this gas preserves its quantum nature at any 
temperature, since only one quantum state is permit- 
ted to be occupied. The number fluctuation of any state 
(the ground state included) is therefore zero. In con- 
trast, the 5hq predicted by the canonical ensemble is a 
smooth non- vanishing function of the temperature, while 
the grand-canonical ensemble still exhibits a divergence. 
This discrepancy between the microcanonical (combina- 
torial) and the other two ensembles remains even in the 
thermodynamic limit. 

One may argue that the Riemann gas is fictitious and 
its spectrum is unrealisable. However, the spectrum 
e n = ln(n) does not increase with n more rapidly than 
n 2 , therefore the existence of a quantum mechanical po- 
tential supporting this spectrum is possible (cf. inverse 
scattering transfor m used in sec tion| lII.B|) . The potential 
has been given in ([Weiss et al\ , \2004) : 



V{x) = Vo In 



(63) 



where Vq and L are positive constants. Within the semi- 
classical approximation the spectrum of this potential is 



V ln(2n + 1) + V In 
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(64) 



where n = 0, 1,. . . and the second term only represents a 
constant energy shift. 



Re cently, LeClair published two works (|LeClairl . 12007 , 
120081 ) developing and applying a finite-temperature field 
theoretical formalism for both boson and fermion gases in 
low spatial dimensions in which he efficiently disentangles 
zero temperature dynamics and quantum statistical sums 
for both the relativistic and non-relativistic cases. His 
alternative approach is base d on an ff-matrix f ormula- 
tion of statistical mechanics (Dashe n et al. , 1969), which 
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redefines the quantum statistical mechanics directly in 
terms of dynamical filling fractions, /(k). Assuming the 
two-body scattering kernel, K, is constant (i.e. constant 
scattering length) he derives, as pedagogical examples, 
the well-known results for the boson 



\C(d/2) 
and also for the fermion gas 
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(65) 



(66) 



In two dimensions the critical temperature for the bo- 
son gas vanishes because of the £(s) divergence at 
s = 1, therefore this dimension needs further con- 
sideration. Due to this instability, LeClair extends 
the examination for energy-dependent two-body kernels, 
K = — •ft(7„/c 2l/ ~ 1 ) (y is a complex number and 7„ is 
constant), for a one-dimensional fermion gas and ex- 
plicitly constructs a quasi-periodic potential, V(x) ~ 
cos (log (x))/x 2cr , in the real space which reproduces the 
given kernel K in the two-body scattering approxima- 
tion. Furthermore, the thermodynamic variables, such 
as density and pressure, are also shown to be physi- 
cally valid (i.e. positive and have finite value) provided 
1/2 < Sft(i/) < 3/2. This fully covers the right hand 
side of the critical strip divided by the critical line, and 
due to the symmetry of £(s) this half-strip can be ex- 
tended to the whole critical strip. His argumentation is 
based on both the non-vanishing, non-divergent nature 
of the physical quantities and also on the assumption 
that an interaction necessarily modifies the thermody- 
namical quantities. If C( y ) would be zero somewhere in 
the critical strip, but off the critical line, then the lead- 
ing order contribution to the thermodynamical quantities 
would not be zero contradicting the original assumption 
- LeClair argues. This contradiction led him to conclude 
that must be non-zero in the 1/2 < v < 3/2 strip, 
which can automatically be extended to the whole critical 
strip by using the symmetries of the Riemann zeta func- 
tion. LeClair, therefore, claims: can have no zeros 
in the given range, consequently the Riemann Hypothesis 
must be true. The basis for this conclusion however, is 
itself an assumption and so does not constitute a proof of 
the Riemann Hypothesis, but does provide another point 
of attack. 

Examination of a similar fictitious, fermionic, many- 
body system has also been considered by Leboeuf and 
lead to the conclusion that "time-periodic dynamical evo- 
lutions have to be considered as serious candidates [f or 
the Hilbert-Polya Hamiltonian]" (jLeboeuf et aLl . l200ll ). 

At the end of this section, let us mention an interest- 
ing interlocking area of statistical physics and number 
theory. A few authors have focused on the connection 
betw een number t heore t ical functions and Browni a n mo- 
tion ( Billingslev . 1973 ; Evangelou and Katsanosl 2005t 



FIG. 9 Function M(n), the cumulative sum of the Mobius 
function is shown with the mean displacement of a random 
walk, ~ \fn for comparison. 
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1998f ) or percolation ( Vardi . 1998). The connection seems 
to be suggestive, especially if one defines the random mo- 
tion through the Mobius function fi(n), i.e. if /i(n) = ±1 
the particle moves up- or downwards, and if n(n) = 
it does not move. Therefore the distance of the particle 
from the origin after n steps is M(n) = X^M^)- The 
importance of this kind of Brownian motion lies in the 
so-called Mertens conjecture. This states if |M(n)| < \fn 
then the Riemann Hypothesis is true. 

Figure [9] shows the path of the particle for the first mil- 
lion steps. Although it is tempting to conclude: the cu- 
mulative sum of n(n) remains bounded by iy/n, this con- 
jecture would actu ally be wrong as te Riele and Odlyzk o 
indirectly proved ( Titchmarsh and Heath-Browni 120031) . 
There is no explicit counterexample known, but we have 
a loose interval [10 14 , ~ 3.6 x 10 10 ] in which there ex- 
ists a n n such that M(n)/y/n > 1 [ Kotnik and Rielel . 
20061 ). Nevertheless, the Mertens conjecture is a suffi- 
ciency condition for the Riemann Hypothesis to be true, 
not a necessary one. Its falsity therefore cannot invali- 
date the Riemann Hypothesis. The failure of the Mertens 
Conjecture at such a high n value, however, does give 
cause for concern regarding numerical evidence for the 
validity of the Riemann Hypothesis. 

However, this is not the only possibility to define a 
random walk either on the £(s) zeros or on the prime 
numbers. In the early 1970's Billingsley defined a ran- 
dom, but finite, walk (|Billingslevl . Il973h based on the 
fundamental theorem of arithmetic. 

Let f(n) denote the number of prime factors of n not 
counting their multiplicity, e.g. /(40) = 2, since 40 = 
2 3 x 5. It can be shown, that on average, numbers below 
N have ln(ln(AT)) factors; a result which on its own is 
a surprise. For example, numbers below e e m jq9566 
have only 10 distinct factors on average. Based on the 
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factorisation one can define the following random walk: 
chose an integer n € [0, N], starting from the origin we 
go up by a unit if 2 divides n and down if it does not, and 
continue the test with 3,4... Although this construction 
does not seem to be as random as a coin-tossing random 
walk and has few flows (e.g. it is biased), Billingsley 
suggests a remedy to these problems and shows how the 
similarity to Brownian motion leads to an Erdos-Kac 
central limit theorem for f(n) 



a< f(n)-Hln(N)) < R " 
. " y/HHN)) - J 
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e- u2 ' 2 du. 
(67) 

Therefore, the probability of f(n) not deviating from the 
expected value ln(ln(iV)) more than a or j3 times the 
standard deviation can be estimated by a Gaussian inte- 
gral. Therefore, the mapping of the number theoretical 
problem onto a Brownian motion helps to derive a limit 
theorem for the number theoretical function f(n). As an 
example, if we chose a = — 1 and /3 = 1 for N = 10 9566 
gives P (-1 < (f(n) - 10) /VW < l) « 0.68, thus ap- 
proximately 70% of the numbers below the chosen N have 
from 6 to 13 distinct prime factors. 



M. Wolf defined random walks in a different way (|Wolf 



19981 ) and could examine the distribution and correla- 
tion of twin-primes (where p and p + 2 are both primes) 
and also of cousin primes (p and p + 4 are both primes). 
He also suggested new random number generators with 
theoretically infinite period based on this kind of ran- 
dom walk, contrary to t he wi dely used random number 



generators ([Press et all 120071) . He also argues and with 



computations demonstrat es the multi -fractal nature of a 
subset of prime numbers ( Wolf , 19891 ). 



IV. CONCLUSION 

'All results of the profoundest mathematical 
investigation must ultimately be expressible 
in the simple form of properties of the integers. ' 

(Leopold Kronecker) 

Since this review is a summary itself in some sense, here 
we only attempt to conclude with some general remarks. 

In many respect the history of the Riemann Hypoth- 
esis is very similar to that of Fermat's Last Theorem, 
which was stat ed in the seventeenth century and solved 
358 years later ( Aczell . 1997 : Ribenboiml . 1993), and along 
the path towards the final proof it inspired and gave birth 
to new areas of mathematics, such as the theory of ellip- 
tic curves. Although the Riemann Hypothesis has not 
been proven or disproven it has already stimulated and 
influenced many areas of mathematics, e.g. L-functions, 
which can be thought of as generalised zeta functions, for 
which a generalised Riemann Hypothesis may hold. In- 
terestingly, for L-functions defined over functional space 



rather than the number field, the similar hypothesis is 
rigorously proven. 

Further evidence also suggests the validity of the Rie- 
mann Hypothesis, let us just think of Levinson's theorem 
guaranteeing that at least one third of the zeros are on 
the critical line. We, however, cannot exclude the pos- 
sibility of the existence of a counterexample to the Rie- 
mann Hypothesis, i.e. a very high lying zero s = a + it 
for which a ^ 1/2. Similarly to the Mertens conjecture, 
the counterexample may occur so high on the critical 
line, that we have no machinery to even calculate zeros 
at that elevation. The immediate impact of such a col- 
lapse of the Riemann Hypothesis would be immense since 
th ere exist numerous "proofs" that are contingent upon 
it ( Titchmarsh and Heath-Brownl 2003 ). 

That said, we cannot miss out in this review one 
computational masterpiece. Not long after World War 
II, in which mechanical and electrical 'computers' were 
often used for encrypting messages (Enigma) and also 
for research (ENIAC), Balthasar van der Pol con- 
structed an electro-mechanical machine which could cal- 
culat e the first few ze ros of the Riemann zeta func- 
tion (jvan der Poll Il947l) . This construction, despite its 
limited achievement, deserves to be treated as a gem 
in the history of the natural sciences. Several decades 
later, on the other end of the spectrum, a state-of- 



by Brent, van de Lune, te Riele and Winter ( 


Brent 


1979: Brent et all 1982; 


van de Lune and te Riele 


. 1983 


van de Lune et al, 1986 


) calculated the first 1.5 


x 10 9 



zeros. Meanwhile Odlyzko (jOdlyzkol 1987 ) explored the 
zeros located around t ~ 10 20 , and showed that all zeros 
(millions of them) he found do exactly lie on the critical 
line. Here we note, that these numerical checking are, of 
their own right, significant achievements, and also have 
influenced the developme nt of fast numerical techniques 
used in physics (see e.g. ( DraghicescuL 19941 Greengardl 
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In this review article we collected a few examples from 
different areas of mathematical physics, starting with 
classical mechanics and finishing with statistical mechan- 
ics, where the Riemann zeta function C( s )j especially its 
zero- and pole-structure, has a highly influential role. 

In the section devoted to classical mechanics, we 
showed how the Riemann Hypothesis can arise in a sim- 
ple mechanical system, a ball bouncing on a rigid wall. 
We also argued how these billiard systems lead to a revo- 
lutionary new way of describing the dynamics of a chaotic 
system by introducing the evolutionary operator. Here 
we also sketched the connection, a trace formula, between 
the dynamics of a chaotic system and the periodic orbits 
of the same system. 

This new descriptive language of dynamics through the 
trace formulae of the Green function is suitable to de- 
velop a new quantisation technique for chaotic quantum 
systems which otherwise was impossible using the stan- 
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dard Bohr quantisation rules. Gutzwiller's trace formula 
has been explicitly mentioned, because the Riemann zeta 
function obeys a very similar expression. Therefore, we 
could compare the two formulae, (j9bf and ([2"5j). and 
imagine what properties a quantum system might have 
if its spectrum mimicked the zeros of the Riemann zeta 
function. 

We also surveyed two other attempts to find a quantum 
system which has a connection to the Riemann zeta func- 
tion. Both of these directions try to associate £(s) with 
the spectrum of the system. The difference between these 
approaches is that one of them relates ((s) to the positive 
energy spectrum, i.e. scattering states, while the other, 
based on the Hilbert-Polya conjecture, proposes systems 
where the negative energies, thus the bound states of 
the system, coincide with the zeros of C( s )- This latter 
case naturally guides us to condensed matter physics and 
statistical mechanics, where one has to evaluate physical 
observable on the lattice points, or derive all thermody- 
namical properties of a given particle-system provided 
the spectrum is given. 

In the sections concentrating on condensed matter 
physics, we first showed how the Riemann zeta function, 
or one of its ancillary functions, arose when we calculated 
the binding energy of a given structure of solid matter. 
Finally we showed how physical requirements for the spe- 
cific heat of a solid can provide zero-free regions for the 
Riemann zeta function. Research in this direction even- 
tually may offer narrower zero-free regions, and comple- 
ment the approach in pure mathematics. 

In the last section, we discussed three main areas of 
statistical physics where the Riemann zeta function and 
its number theoretical aspects influence the behaviour 
of a physical system. Firstly, we considered the low- 
temperature phase transition of bosons, and showed that 
the pole structure of £(s) prohibits Bose-Einstein con- 
densation in one- and two-dimensional uniform systems. 
We also reviewed the 'grand canonical catastrophe' of 
an ideal Bose gas, where the predictions of two ensem- 
bles widely used in statistical physics contradict each 
other, showing, therefore, that these ensembles cannot 
be equivalent to one another. Finally, we examined a 
possible Brownian motion model for the number theoret- 
ical Mobius function. 

It would not be without precedent if a completely new 
theory or a new mathematical language is needed in 
which the Riemann Hypothesis can be 'worded' naturally 
for the hypothesis to be finally proved. As has happened 
earlier with mathematics, natural science, and in particu- 
lar physics, can give impetus and motivate new directions 
perhaps leading to the final proof. It is amazing and cap- 
tivating to see that a purely number theoretical function 
has so many direct links to classical and modern physics. 

Nowadays we are not surprised by Galileo's famous 
keynote: "[Nature] is written in the language of mathe- 
matics, and its characters are triangles, circles, and other 



geometric figures without which it is humanly impossible 
to understand a single word of it; without thes e, one wan- 
ders about in a dark labyrinth" ( Drakel Il957t ). Probably 
we are not meandering in a labyrinth, but we are defi- 
nitely puzzled by the overwhelming difficulty of proving 
the Riemann Hypothesis. We simply do not know as 
yet whether physics will ultimately help in understand- 
ing such an elegant mathematical statement as the Rie- 
mann Hypothesis, but we are definitely witnessing the 
intertwining and invigoration of both disciplines. The 
authors can only express their hope that this work has 
to some extent captured the imagination of the Reader 
and. if so, it has fulfilled its intended aim. 
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